To appear in Phys. Rev. D 



Preprint, March 1995, 



gr-qc/9503061 



Exact Hypersurface-Homogeneous Solutions in Cosmology and Astrophysics 

Claes UgglaB 

Department of Physics, University of Stockholm, Box 6730, S-113 85 Stockholm, Sweden 

Robert T. Jantzen 

Department of Mathematical Physics, Villanova University, Villanova, PA 19085, USA and 
International Center for Relativistic Astrophysics, University of Rome, 1-00185 Roma, Raly 



in 



KjeU Rosquistm 

Department of Physics, University of Stockholm, Box 6730, S-113 85 Stockholm, Sweden and 
RGGR group, Chimie-Physique, Universite Libre de Bruxelles, Brussels, Belgium 



o 



A framework is introduced which explains the existence and similarities of most exact solutions of the Einstein equations with 
a wide range of sources for the class of hypersurface-homogeneous spacetimes which admit a Hamiltonian formulation. This 
class includes the spatially homogeneous cosmological models and the astrophysically interesting static spherically symmetric 
models as well as the stationary cylindrically symmetric models. The framework involves methods for finding and exploiting 
hidden symmetries and invariant submanifolds of the Hamiltonian formulation of the field equations. It unifies, simplifies and 
extends most known work on hypersurface-homogeneous exact solutions. It is shown that the same framework is also relevant 
to gravitational theories with a similar structure, like Brans-Dicke or higher-dimensional theories. 
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I. INTRODUCTION 

Exact solutions have always played a central role in 
the investigation of physical theories whose content is en- 
coded in a complicated set of differential equations. The 
theory of general relativity is indeed an example of this. 
A number of exact solutions of Einstein's field equations 
have been of key importance in the discussion of physi- 
cal problems. Solutions have been found which describe 
black holes, stellar interiors, gravitational waves, and 
even the large scale structure of the universe itself. Exact 
solutions have also served as a guide to point out mathe- 
matical features of the theory The Taub-NUT-M (Taub- 
Newman-Unti-Tamburino-Misner solution, for ex- 

ample, has been of crucial importance for the very def- 
initions one uses in describing the singularities of the 
full theory. Thus exact solutions may point out features 
which are not just special to themselves but character- 
ize in some way properties of a wider class of solutions. 
They may also play a role as "building blocks" for more 
general solutions. For example, in certain ways the gen- 
eral spatially homogeneous cosmological model near an 
initial singularity can be understood in terms of very spe- 
cial exact solutions, notably the Kasner and the vacuum 
Bianchi type II solutions, which to some extent also de- 
scribe aspects of general cosmological singularities 
Sometimes exactly solvable problems are even used as a 
guide in developing ideas for the construction of more 
general theories, e.g., quantum gravity. For example, 
solvable problems in spatially homogeneous (SH) cosmol- 
ogy have been used to implement a number of different 
quantization schemes. 



Thus there is ample motivation to try to find exact so- 
lutions. Indeed, the book by Kramer et al Q is largely 
dedicated to the listing of exact solutions. Several chap- 
ters of that book deal with hypersurface-homogeneous 
(HH) solutions, a class for which the Einstein equations 
reduce to more manageable ordinary differential equa- 
tions. Within the class of HH solutions there are several 
subclasses of considerable physical interest, the cosmo- 
logical SH models and the astrophysical static spherically 
symmetric spacetimes being the most prominent ones. 

Since the birth of general relativity nearly 80 years 
ago, an overwhelming number of HH solutions have been 
produced. A look at physics abstracts shows that this 
production continues even today at a considerable pace. 
However, often this search is undertaken as an end in 
itself without attempting to understand how particular 
successes fit into a larger scheme and without employ- 
ing any systematic method of attack revealing possible 
underlying mechanisms. Exceptions do exist though, as 
illustrated very nicely by the numerous approaches to the 
problem of finding vacuum solutions for spacetimes with 
one or two commuting Killing vector fields [^jj^. How- 
ever, these techniques have not contributed much to the 
problem of finding HH solutions. 

It is the purpose of this article 

• to give a basic underlying explanation of why ex- 
act solutions arise for models with nonnuU homo- 
geneous hypersurfaces, 

• to provide a set of techniques which makes it pos- 
sible to obtain these solutions in as simple a form 
as possible, 
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• to show how these techniques are apphcable to a 
wide set of different physical problems, and 

• to unify an otherwise seemingly unrelated zoo of 
particular results. 

Rather than exhaustively treating all possible 
wide variety of examples will be discussed which illustrate 
the utility of the approach presented here. The choice 
of examples is a subjective one reflecting our particular 
tastes. Among the many models and source types con- 
sidered, the exact SH perfect fluid solutions will be more 
exhaustively surveyed, updating the work of Kramer et 
al. Scalar fields, not considered in that catalog but now 
currently fashionable and physically interesting, will also 
be examined. Special attention will also be given to the 
static spherically symmetric perfect fluid models which 
are important in astrophysical applications. 

The present approach will use a Hamiltonian/Lagran- 
gian formulation of the field equations, using the Hamil- 
tonian constraint together with the Lagrangian equa- 
tions. This may come as a surprise considering the state- 
ment by Kramer et al (p. 131) that this formulation 
"is not well-adapted to searching for exact solutions. . ." , 
but upon second thought, this is quite natural since the 
Hamiltonian function contains all the dynamical content 
of the Einstein equations. This enables one to study a 
single function instead of a whole system of equations, 
armed with many powerful techniques from classical me- 
chanics. The economy of the Hamiltonian approach also 
reveals the close mathematical relationships among dif- 
ferent types of models. These relationships are often 
obscured by the particular way in which a particular 
physical problem suggests expressing the field equations. 
For example, the usual approaches to static spherically 
symmetric models and SH cosmological models are quite 
different, but the Hamiltonian formulation shows their 
mathematical similarities. Furthermore, the Hamilto- 
nian approach used in this article will also show how 
most exact solutions are either associated with essentially 
1-dimensional problems (in terms of degrees of freedom) 
or with problems which admit a sufficient number of a 
certain type of Hamiltonian symmetry. 

The article will proceed as follows. Since the Hamilto- 
nian is essential for this discussion, Hamiltonians will be 
derived for a variety of different problems in section 2. 
This section starts with an outline of the way in which 
the Hamiltonian function is obtained for the HH mod- 
els. This is followed by the explicit evaluation of this 
function, considering separately the diagonal and nondi- 
agonal models. For the more numerous diagonal models, 
the gravitational Hamiltonian is discussed first and then 
a number of source contributions to the Hamiltonian are 
examined. For the fewer nondiagonal examples, each case 
is studied individually, and only the vacuum and perfect 
fiuid Hamiltonians are derived. In all cases the Hamilto- 
nian is put into a "canonical" form that makes evident 
the mathematical similarities among different HH mod- 
els. 



In section 3 the generalized Friedmann equation is re- 
viewed. This equation often arises in the context of 1- 
dimensional invariant submanifolds and as a part of solv- 
able higher-dimensional problems. It therefore frequently 
plays an important role as a fundamental building block 
when it comes to finding exact solutions. 

To solve a problem with two or more degrees of freedom 
it is necessary to find a sufficient number of symmetries 
which makes it possible to decompose the problem into 
smaller solvable parts. In section 4 a particular type of 
Hamiltonian symmetry is discussed, a so-called Killing 
tensor symmetry. It turns out that this symmetry, which 
usually is "hidden," is responsible for all known solv- 
able Hamiltonian problems with two or more degrees of 
freedom that we are aware of. This symmetry is a gener- 
alization of symmetries related to the existence of cyclic 
variables and Hamilton- Jacobi separability. To exploit 
the Killing tensor symmetries in order to find explicit 
solutions, one has to find symmetry compatible depen- 
dent and independent variables. A way of doing this is 
presented. 

Section 5 presents one method for finding invariant 
Hamiltonian submanifolds, i.e., consistent subsystems of 
the Hamiltonian equations. This in turn leads to partic- 
ular solutions. The existence of invariant submanifolds 
is an important issue in the search for exact solutions 
since it is usually impossible to solve the most general 
problems. 

Section 6 lists problems leading to exact solutions and 
indicates how to solve them by referencing the relevant 
parts of this article. The models in this section are again 
divided into diagonal and nondiagonal models. However, 
the diagonal models previously arranged in section 2 ac- 
cording to how their Hamiltonians are evaluated are now 
instead collected together according to the dimension of 
their intrinsic symmetry group, i.e., the group of sym- 
metries of the intrinsic geometry of the individual ho- 
mogeneous hypersurfaces. This categorization reveals an 
underlying mathematical unity of entire classes of phys- 
ically distinct models and allows them to be treated col- 
lectively. 

The methods developed in this article are also appli- 
cable in contexts beyond 4-dimensional classical general 
relativity. This is discussed in the concluding section 
7. As examples some remarks are made about higher- 
dimensional theories, nonminimally coupled scalar field 
theories, and quantum cosmology. 

In section 8 the present approach is compared with 
some other exact solution techniques, followed by a con- 
cluding discussion which addresses a variety of other is- 
sues. 



II. A HAMILTONIAN APPROACH TO THE 
FIELD EQUATIONS 

In this section the Hamiltonian will be evaluated for a 
wide variety of HH models. However, the goal is not just 
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to produce a Hamiltonian, but to obtain it in a certain 
"canonical" form which reveals the mathematical simi- 
larities between different HH models. This is not only 
useful for the present purpose of finding exact solutions, 
but may also serve as a starting point for qualitative 
analysis of the many problems which cannot be solved 
exactly. Furthermore, one can study a single function 
rather than be overwhelmed by a whole system of equa- 
tions which also hides the mathematical connections be- 
tween different kinds of models. This section concludes 
with a discussion of the general form of the Hamiltonians 
which have been obtained in the individual cases. 



A. The line element 

To explicitly obtain a Hamiltonian for a given model, 
it is necessary to introduce a line element. In this article 
we will consider the Taub-NUT-M model, to be discussed 
below, and certain models with nonnuU homogeneous hy- 
persurfaces. For the latter models the line element can 
be expressed in the form 

ds^ = eN{XfdX^ + gab{X)oj''uj'' , (2.1) 

where the 1-forms (a = 1,2,3) are associated with 
the symmetry group acting on the homogeneous hyper- 
surfaces. The quantity e = UaTi" is the sign of the norm 
of the unit normal to the homogeneous hypersurfaces, 
having the value —1 for the SH models and 1 for the 
static ones (a assumes the values 0,1,2,3 with refer- 
ring to the component associated with the A direction). 
For the SH models the single independent variable A is 
a time variable t and N is the familiar lapse function. 
In the static case A is instead a spatial variable. The 
gravitational degrees of freedom are associated with the 
component functions gab- 

For later purposes it is convenient to introduce the 
function x defined by 

X - A^(taub)/A^ , %aub) = 12|ff|'/' , (2.2) 

where g — det(gah) = The function x is the re- 

ciprocal of the relative "slicing gauge function" Af — 
-^/-^(taub) which is sometimes used in 3 -t- 1 canonical 
gravity as a relative lapse. One must fix or a; to deter- 
mine the parametrization of the family of homogeneous 
hypersurfaces by the independent variable A; this will be 
referred to as a choice of slicing gauge. The notation 
-^(taub) comes from the very useful slicing gauge intro- 
duced by Taub ^ in the context of SH cosmology. 

The HH models admit simply transitive or multiply 
transitive (MT) homogeneity groups. The models which 
admit a simply transitive three-dimensional homogeneity 
group are called Bianchi models, for which the 1-forms 
may be chosen to satisfy 

dw" = -iC%cw'' Aw", (2.3) 



where C^tc are the components of the structure constant 
tensor of the Lie algebra of the homogeneity group. The 
Bianchi models are divided into two classes, class A and 
class B, according to the vanishing or nonvanishing of the 
trace C^ab- For the class A models one can choose C^bc = 
'n^°''^€abc, where the parameters n*-"^ characterizing the 
various symmetry types can be chosen to have the values 
in Table 1; explicit coordinate representations of the 1- 
forms can be found in [Q. Class B models which admit 
a Hamiltonian formulation will be discussed below. 

Among the Bianchi models there is a special family of 
Bianchi type I, II, III, V, VHq, VII,,, VIII, and IX models 
which admits MT symmetry groups. There are also MT 
models which do not admit a simply transitive subgroup 
(i.e., they are not Bianchi models); the most prominent 
ones are the SH Kantowski-Sachs models and the static 
spherically symmetric models. 

B. How to find the Hamiltonian function 

Not all HH models admit a Hamiltonian formulation 
of the field equations. However, of those that do, the 
Lagrangian and Hamiltonian have the usual simple forms 
in terms of kinetic and potential terms 

L = T-U, H = T + U. (2.4) 

The Hamiltonian must vanish as a consequence of one of 
the field equations; this is the "Hamiltonian constraint" 
H = 0. From the expression for H in terms of the veloc- 
ities, one can read off the Lagrangian as well by a simple 
change of sign. The Hamiltonian constraint and the La- 
grangian equations are both needed since a mixture of 
first and second order differential equations is often re- 
quired to find explicit solutions. Since the Hamiltonian 
constraint is of crucial importance in obtaining solutions 
and is the most economical way of describing the field 
equations, our method will be referred to as a Hamilto- 
nian approach. 

The starting point for obtaining the above Hamiltonian 
is the ADM Hamiltonian of the 3 -I- 1 or ADM approach 
described in Misner, Thorne and Wheeler ||]. This 
approach reformulates the Einstein equations expressed 
with respect to a spacelike slicing as a parametrized 
Hamiltonian system with constraints. Since all their for- 
mulas involve explicitly the sign e of the norm of the unit 
normal to the slicing, they continue to hold for a timelike 
slicing. The ADM Hamiltonian is a linear combination of 
the super-Hamiltonian and supermomcntum constraint 
functions (whose vanishing is equivalent to the 00 and Oa 
components of the Einstein equations respectively) , with 
the lapse and shift as the respective Lagrange multiplier 
coefficients. However, in the present application the su- 
permomcntum constraints can be solved and the shift 
set to zero, reducing the ADM Hamiltonian to the lapse 
times the super-Hamiltonian. Allowing the lapse func- 
tion in this Hamiltonian to depend explicitly on the grav- 
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itational and possible source variables leads to Hamilto- 
nian equations which differ from each other by terms in- 
volving the super-Hamiltonian, required to vanish as the 
remaining unsolved constraint. 

This specialized ADM Hamiltonian is the first step to- 
wards obtaining the desired Hamiltonian for HH models 
and will be referred to simply as the Hamiltonian func- 
tion. It has the following form |8| 



H, 



(G) 



(source) ; 



(2.5) 



where Gap is the Einstein tensor and T^/j is the total 
energy-momentum tensor, including a possible cosmolog- 
ical constant term, while k is the appropriate multiple of 
the gravitational constant. The field equations require 
H = Q. 

The gravitational part of H is 



^^(G) - ^(G) 



N\g\^'^[K\K\-K\K' 



:(3) 



R] 



= 12x-' \g\[K\K'' a - K\K\ + e '^^^R] , (2.6) 

where the extrinsic curvature tensor has the usual ex- 
pression 



-\N-^g--g,, 



-kAaV'^rac (2.7) 



(a dot indicates the A derivative) and is the scalar 
curvature of the homogeneous hypersurfaces. 

To go further one must discuss how the source en- 
ters the Hamiltonian. Rather than attacking the general 
case, we will confine our attention either to a scalar field 
which contributes one additional degree of freedom, or to 
sources whose additional degrees of freedom can be elim- 
inated by integrating the source field equations in such a 
way that the source variables can be expressed entirely 
in terms of the gravitational variables and constants of 
integration. The second step towards obtaining the HH 
Hamiltonian involves inserting those expressions for the 
source variables into the Hamiltonian function. 

The third step is to solve the supermomentum con- 
straints, if nontrivial. If these latter constraints are holo- 
nomic, i.e., can be expressed as the vanishing of a total 
derivative, they can be integrated. In this case they can 
be used to reduce the number of gravitational degrees of 
freedom, and the above Hamiltonian function gives the 
correct field equations. On the other hand if these con- 
straints are nonholonomic, the reduction of the number of 
degrees of freedom may not be possible or may not yield 
a Hamiltonian which gives the correct field equations be- 
cause of symmetry considerations Q. Because of these 
difficulties one should always check that the Hamiltonian 
function yields the correct field equations. 

Recently the relationship between topology and the 
Hamiltonian formulation of the field equations has at- 



tracted some interest | p^ , pl| . If the symmetry type al- 
lows a closed topology for the homogeneous hypersur- 
faces, then a Hamiltonian formulation exists. However, 
the nonexistence of such a topology does not exclude a 
Hamiltonian formulation, as exemplified by certain class 
B type VI?i models treated below. 

However, one does not just want to obtain a Hamilto- 
nian for the problem, but also to express it as simply as 
possible. This is helpful not only in finding exact solu- 
tions but also for a qualitative analysis of dynamics which 
cannot be described by exact solutions. The fourth and 
final step towards obtaining the HH Hamiltonians in a 
useful form is accomplished by choosing new variables 
which diagonalize the kinetic part T of the Hamiltonian. 
The way in which this is achieved differs for models which 
are diagonal or nondiagonal, i.e., models for which the 
line element can or cannot be expressed in diagonal form. 
These two cases will be treated separately. 



C. Diagonal models 

The line element for diagonal HH models can be writ- 
ten in the form 

ds" = eN{\f d\^ + .9ii(A) [u^f + (?22(A) {u?f 

+.g33(A) [J'f . (2.8) 

It is convenient to introduce the following variables 



(2.9) 

The diagonal models all have a Hamiltonian formula- 
tion. The kinetic part T{q,) of the gravitational Hamilto- 
nian can not only be diagonalized but even made "confor- 
mally flat" by a linear transformation of the /J" variables. 



1. Class A models 

For these models the momentum constraints are satis- 
fied identically. The parametrization 




■ /30 + /3+ + V3/3- 
13" + 13+ - V3/3- 
13° - 2(3+ 



(2.10) 



introduced by A^Iisncr ||l^ , leads to the "conformally flat" 
form 



(2.11) 



of the expression (2.6) for the kinetic part of the gravi- 



tational Hamiltonian. Here ijab = V (^, -B = 0, -f , — ) 
are the orthonormal components of the Minkowski metric 

-VOQ = V++ =V — = 1 (2-12) 
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and the quantity x is the conformal factor, mak ing it 
more convenient to use than N . Equation ( ^.6| ) also 
yields the gravitational potential 

a=l 

-2n(2)„(3)g22g33 - 2^(3)^(1)533511] • (2.13) 

In addition to being mathematically convenient, the 
variables are also adapted to physical quantities. The 
quantity \g\^/'^ — — e^^ is the volume-element 

density. The A derivatives of the (3'-', (3^ variables are re- 
spectively proportional to the expansion and the nonzero 
shear components of the congruence normal to the homo- 
geneous hypersurfaces. 

The possibility of choosing the timelike direction along 
different axes in the static case yields an abundance of 
models. To avoid being lost in details, the signature free- 
dom will be restricted by requiring lo^ and cj^ to be space- 
like, a condition which will be assumed in the rest of this 
article. This leads to the line element 



(2.14) 



With this choice of timelike direction, the gravitational 
Hamiltonian takes the form 



i?(G) - ixr;^B/3^/3^ + C/(G) 

= ^xr^ABP^P"" + 12x-^e^^"v*i(3^) , (2.15) 



where 



V* 



(2.16) 



An interesting subclass of these models is obtained by 
setting /3~ = (equivalent to i?i = i?2). For the Bianchi 
types I, II, VIIo, VIII, and IX, these models admit multi- 
ply transitive symmetry groups which correspond to local 
rotatio nal sy mmetry for the choice made for the timelike 
axis in (2.14). The Bianchi type I and VIIq models of this 
type coincide and are plane-symmetric. Only the type I 
models will be referred to in what follows. 



2. Bianchi type V and Vlh models 

Most class B models do not admit a Hamiltonian for- 
mulation, but some particular nondiagonal and all di- 
agonal models (consistent with the field equations and 
sources studied here) do. In the following we will only 
consider SH class B models since the static class B mod- 
els seem to lack physical interest. The SH class B Bianchi 



types V and Ylh admit diagonal cases provided that the 
structure constants have the form iTUlSjlitl 



32 



-hq^ ■ (2.17) 



The line element has the same form ( 2.14 ) as in the class 
A case, but the variables Ra are related by an algebraic 
constraint obtained from the single nontrivial holonomic 
momentum constraint. Canonical choices for the struc- 
ture constants are 



V 
VI^ 



9 = 1, 



a = 1 , 
a > 0. 



(2.18) 



Explicit coordinate representations of the 1-forms may 
be found in @,|l|. 

The following parametrization solves the momentum 
constraint and diagonalizes T(g) 



/ (3° - c{q-3a)f3' 
.13^/ \ ^o + 2cg/3>< 



(2.19) 



where c = {q"^ + 3a^)~^^^. Note that for the Bianchi type 
V case, this reduces to the class A parametrization with 
f3^ = and (3~ = (3^ . The parameter combination cq 
takes values in the interval [0, 1) with the endpoint value 
1 corresponding to the class A limit of type VIq which 
arises from a — > 0, with q ^ Q. The parametrization then 
reduces to the corresponding class A case with j3^ ~ 
and /?+ = — /3^. It is sometimes convenient to choose 
a = 1 in Bianchi type V, leading to = 3, since this 
choice is the one usually used for the open isotropic FRW 
models which are obtained by setting = in type V. 

With the above parametrization the gravitational 
Hamiltonian has the expression 



(G) 



(2.20) 



3. Multiply transitive models 

The MT models we have encountered so far are those 
belonging to class A and the isotropic open type V model. 
There are some additional ones described by the same 
form (2.14) of the line element as in the class A case with 
the Misner parametrization restricted by the condition 
/3~ = corresponding to i?i = i?2- For this case, the 
{lo^ Y' + {uj'^ Y PEH't of the line element represents a 2-space 



of constant curvature a € {1,0,-1}, while 



,3 - 



dx^ is 



an exact differential. Explicit coordinate representations 
of the other two 1-forms may be found in 

For this class of models the case a — coincides with 
the locally rotationally symmetric (LRS) Bianchi type 
I models already considered. The case cr = — 1 is an 
LRS class B Bianchi type HI case. The case cr = 1 is 
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the exceptional case where no three-dimensional simply 
transitive subgroup exists, not falling under the Bianchi 
classification. The SH models of this latter case are called 
Kantowski-Sachs (KS) models, while the static models 
are the spherically symmetric models. 

The 3-curvature '■^•'i? reduces to the curvature of the 
constant curvature 2-spaces '^^^R ~ Ri^^cr = e^"^^ a, 
leading to the gravitational Hamiltonian 

a;(_^02 + ^+2) _^ 24eax-'e*f"-^f'^ . (2.21) 



(G) 



The type III models can be obtained from the previ- 
ously discussed type VI/i models. The values a — 1 = q, 
for which c — ^, correspond to the usual Bianchi type III 
= VIft,=_i structure constants. However, the equivalent 
choice a — = q leads to c = 1 and the same Hamilto- 
nian as (2.21) (with tr = — 1) for this case. In both cases 
one has cq = ^. The different values of c are related by a 
translation in the (3^ variables associated with a scaling 
of the corresponding 1-forms. 

When not otherwise stated a reference to MT mod- 
els will refer to the LRS type I, HI models and the KS 
and static spherically symmetric models described by the 
Hamiltonian (2.21). 



4- Source terms 

Any combination of potential terms of the following 
form may be considered, representing different arrange- 
ments of sources. 

a. Cosmological Constant. A cosmological constant 
term in Einstein's equations will be considered as an ad- 
ditional term {—A/ K)ga/3 in the total e nerg y- momentum 
tensor Tq,/?- Inserting this term in Eq. ( |2.5|) gives rise to 
the potential 



U, 



(A) 



(2.22) 



b. Scalar Field. A scalar field adds another depen- 
dent variable and contributes both to the kinetic part of 
the Hamiltonian and to the potential |l3] 



(sc) 



ix/3t2, = -2te-ie6^V(,e)(/3^), (2.23) 

where V(sc)(/3^) is the scalar field potential and where 
= y/ K/6(f) is the relation between the scaled scalar field 
variable and the usual scalar field (j). It will be con- 
venient to introduce = (3^ and an exponential scalar 
field variable analogous to the gravitational variables Ra 
by defining 



-R4 



(2.24) 



When convenient the index a will take the values in the 
set {1,2,3,4} and the range of the index A may include 
the value "f." 



c. Perfect Fluids. For a perfect fluid source with en- 
ergy density p, pressure p, and 4- velocity u", the energy- 
momentum tensor is 



(2.25) 



For the cosmological case e = — 1, the 4- velocity for a 
diagonal source must be u" — n". For the static case 
e = 1, the 4- velocity must be alo ng th e timelike third 
direction Ua = u^S^a (recall Eq. ( 2.14 )). According to 
Eq. (2.5), the fluid contributes the potential 



2KN\g\^^^p (SH case) 
2KN\g\^^'^p (static case) 



(2.26) 



to the Hamiltonian but an equation of state is needed to 
express this entirely in terms of the gravitational vari- 
ables using the conservation equations. 

In the SH case an equation of state p — p{p) is imposed. 
Note that dust models do not exist in the static diagonal 
case without additional source terms (to prevent the dust 
from collapsing) . Since the fluid potential is proportional 
to p in this case, an equation of state p — p{p) is therefore 
imposed instead. 

Following Misner, Thorne and Wheeler it is conve- 
nient to introduce the baryon number density n and the 
chemical potential p = {p + p)/n which satisfy 

d\nn= {p + py^dp, d\n p = {p + p^^dp . (2.27) 

The conservation equations then imply 

n"Tj,fj = ip + p){u^[u^,[3 + ilnnp),f3U^] 

+Ua;/3u'' + (lnAi);a}n" = 0. (2.28) 

In the SH case the conservation equation reduces to 







ng 



1/2 



const . 



(2.29) 



This equation together with ( ^.27 ) leads to p p(/3") 
which gives the fluid potential 



(fluid) 



= 2AKx-'e'^''piP°) . 



(2.30) 



For the static case the conservation equation reduces 
to the simple relation 



dln/i = 



dp 



P + P 

which may be integrated to yield 



pRs — pe'^ — const , 



p (X e 



-0' 



which means that p = p{0^) = p(/?° — 2/3^ 
the the potential 

?7(fl,id) = 24Kx-ie*^'^°p(/3° - 2/3+) . 



(2.31) 

(2.32) 
This yields 

(2.33) 
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For a fluid with the equation of state p — ( 7 — 
implying /i — ^p/n, the above relations (2.27) can be 
integrated to yield 



P/P(0) = {n/n(a)y , = Wrt(o)) 



7-1 



(2.34) 



and without loss of generality one can set p(o) — "-(o) 
leading to 



(2.35) 



The values 7 = 1, |,2 describe respectively dust, radia- 
tion, and stiff perfect fluids. 
For the SH case, this leads to 

P = Pio)e''''"'\ (2.36) 
where the choice pjo) = P leads to the fluid contribution 



(2.37) 



to the Hamiltonian. 

In the static case one often assumes p — p(Q) + {rj—l)p 
as an equation of state, where 77 and /9(o) ^-re constants. 
This includes both p = (7 — l)p when p(o) = 0, so that 
V = 7/(7 ^ 1) J &s well as the case of constant energy 
density obtained by setting 77 = 1. Inserting this into 
Eq. (|2^3l|) yields 



iP{0) +P(0)) 



V P{o) , 



(2.38) 



where p(o) is a- constant of integration. This leads to the 
fluid potential 



P{0) 



=6/3" 1 



(2.39) 



Occasionally, particularly in the SH case, one considers 
several noninteracting perfect fluids with a common 4- 
velocity. In this situation one has a perfect fluid potential 
for each such component fluid. 

d. Electromagnetic Fields. An electromagnetic field 
has the energy-momentum tensor B 



(2.40) 



The simplest electromagnetic fields are aligned with one 
of the axes, leading to a diagonal energy-momentum ten- 
sor. Examining Maxwell's equations for an electric or 
magnetic field or some combination of the two along a 
single axis, one finds which directions are allowed (if any). 
For the class A and MT models the third axis is such a 
direction, leading to an energy-momentum tensor with 
nonzero components 



T% = T^ = -T\ = -T\ = -p. 



(2.41) 



An examination of the conservation equations similar to 
that of the perfect fluid case shows that the positive quan- 
tity 



911922P - 



-1 2 
e 



(2.42) 



is a constant. The electromagnetic contribution to the 
Hamiltonian then takes the simple form 



U 



(2.43) 



In the static cylindrically symmetric case the inho- 
mogeneity together with the singular axis of symmetry 
breaks the duality symmetry, and a magnetic field must 
lie along one of the spatial directions not associated with 
the spatial independent variable A. However, the poten- 
tial is identical with the previous case once a suitable 
permutation of the dependent variables Ra is made. 



D. Some nondiagonal models 

Similar methods can be used to treat some of the non- 
diagonalizable cases. Some of the simplest nondiagonal 
cases with a Hamiltonian formulation will be considered 
to illustrate this procedure. Solving the momentum con- 
straints leads to a problem with an additional nondiag- 
onal gravitational degree of freedom which is associated 
with a cyclic variable. Using the associated constant mo- 
mentum, one obtains a reduced Hamiltonian for the di- 
agonal degrees of freedom with an effective potential left 
behind in the kinetic part of the Hamiltonian. This is 
analogous to the centrifugal potential which appears in 
the central force problem. 



1. The Taub-NUT-M model 

The Taub-NUT-M model is an excellent example to 
study in this context for two reasons. First it has a ho- 
mogeneous slicing whose causal character changes from 
spacelike to timelike and back to spacelike again. Second 
it is a nondiagonal model which pieces together diago- 
nal models of each causality type, illustrating the way in 
which some nondiagonalizable models behave. 

Its line element can be put in the form 



(2.44) 



where the are the same 1-forms as for the diagonal 
type IX models and z, g, and w are functions of A. The 
function z is a slicing gauge function. The function g 
is positive in the SH Taub region, negative in the static 
NUT region, and zero at the bridging null hypersurfaces 
between them. The A coordinate lines are null in contrast 
with the usual orthogonal coordinate lines for which the 
line element is diagonalized in the Taub and NUT re- 
gions. 

Due to the existence of null hypersurfaces, it is per- 
haps most straightforward to specialize the full curvature 
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scalar Lagrangian to this case, removing a total A deriva- 
tive to obtain a Lagrangian vahd for the entire spacetime. 
This is equivalent to the ADM Lagrangian in the sep- 
arate Taub and NUT regions with the shift and lapse 
freedom fixed by the null condition on the A coordinate 
lines, modulo a conformal factor z representing the free- 
dom remaining in the parametrization of the slicing. The 
corresponding Hamiltonian is given by 



2) = 0. 



(2.45) 



The kinetic part can be easily diagonalized for 5 7^ and 
(7 = separately, but not for all values of g simultane- 
ously which is needed to describe the full Taub-NUT-M 
spacetime. 



2. Stationary cylindrically symmetric models 

The stationary cylindrically symmetric models have a 
line element which can be written as Q 

ds^ = N^dX"^ - e^^\dt + Cd(j)f + e^'^'d^^ + e^'^'dz^ , 

(2.46) 

where /3", N and C are functions of the independent 
variable A which is interpreted here as a radial coordinate 
ordinarily denoted by the symbol p. As for the diagonal 
cases, the Hamiltonian is given by 

H - -2iV|g|i/2^"n'3(G„^ - kT^^) = . (2.47) 

Expressing the variables fi"^ in terms of the Misner 
parametrization, the vacuum Hamiltonian assumes the 
explicit form 

H = ix(77AB/3^/3'' + T^e^^"' C^) = . (2.48) 

The momentum pc associated with the cyclic variable C 
is constant, leading to 

C = Ue'-'^'^'pc (2.49) 
and the reduced Hamiltonian 

H = ^xt^abP^P^ + 12x-^ e-'^^'^' pc^ = . (2.50) 



C^bc = n'^°'^eabc + cl5'^z (no sum on a) , 

„(i) = -n(2) = l, n(3)=0, a=\. (2.52) 

The 3-metric gab can be conveniently parametrized by 
gab = S\S\g'^^, where g'^^, = diag(e2'3\ e^/^', e^/^') and 
the matrix S'^f, lies in the 3-dimensional special automor- 
phism group of the type VI_i/9 Lie algebra 



cosh6l3 -sinh^^ \/26l2 
(S"'b)=( -sinh6'3 cosh6'3 ^/2e^ 
1 



(2.53) 



The momentum constraints remain to be satisfied 
in this nondiagonal case. Using the standard Misner 
parametrization for the /3" variables and the above 0" 
variables in the constraint equations leads to |l^ 

6»2 = 61^ = -13+ . (2.54) 

These conditions together with the transformation 

leads to the Hamiltonian 



13'^ = P\ 13+^-^(3+^ e^=^e^ 



(2.55) 



[32e4^""2v^'3+ ^ 2APe^^^~''^^"] , (2.56) 

where £ is the same constant as in the diagonal case . 
Since the fluid is orthogonal, it follows from Eq. ( |2.5| ) 
that it has the same potential as in the diagonal case. 

The momentum p^ associated with the cyclic coordi- 
nate if is constant, leading to 



ip = Ix-'e'^^^p^ 



(2.57) 



and the reduced Hamiltonian 

H = ix[-(/30)2 -I- (/3+)2] + ^-i[|e4V3^+p^2 

^32g4^°-2v^r + 24re3(2-7)/3°] = . (2.58) 

Note that ip — = p^ reduces this case to the corre- 
sponding diagonal case. 



4- Spatially homogeneous class A models belonging to the 
symmetric case 



3. Spatially homogeneous Bianchi type V7_j/g models 

The orthogonal perfect fluid models (u" = n") of this 
type permit a nondiagonalizablc line element of the form 

ds^ = -N{tfdt^ + 5afc(<)w"w'' , (2.51) 

where are the invariant 1-forms for this symmetry 
type, with structure constants 



The class A perfect fluid models with an equation of 
state p = (7 — l)p admit the special case where the fluid 
4- velocity has a single nonzero spatial component W3 and 
where gab in the line element 

ds^ = -N{tfdt^ + gab{t)uj''i^^ (2.59) 

has only one nonzero offdiagonal coefficient 512 — 521- 
This is referred to as a "symmetric case" . For Bianchi 
type II it is now more convenient to choose the alternative 
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structure constant values {n^^\n^^\n^^'>) — (1,0,0) in 
place of those of Table 1. 

As in the previous case, gat can be conveniently 
parametrized by gab = S'^aS'^bg'cdi where is the same 
but the special automorphism matrix S'^b is now 

/ C3 -nWs3 0\ 
iS\) - n(2)53 C3 , (2.60) 
V 1/ 

where 

e"' = 2-i/2[(n(i))2 + (n(2))Y/^ 

C3 = coshm(3)6'3 , S3 = (m^^))"^/^ sinhm^^^ei^ . (2.61) 

For Bianchi types I and II where some of these expres- 
sions are undefined, one defines them as the limit in which 
n^^^ and in Bianchi type I one then lets n*^^^ ^ [Q. 
The Hamiltonian constraint (|]^) is then 

H = r(G) + t/(G) + [/(fluid) = , (2.62) 
where the kinetic part is given by |Q 

= HvabP^P'' + ^^e-^^\h^f{ey] , (2.63) 

while the gravitational potential C/(g) is the same as for 
the diagonal case. The fluid potential is given by 

[/(fluid) = 2At/Vgi/2n"r/T„0 

= 2KNg^'\p + p){n"u^f -p] 
^2'iKx-'^e^^''[{p + p)Y ~p], (2.64) 

where Y = {n°'Ua)'^ . As for the diagonal case for an 
equation of state p ^ (7 — l)p one can introduce p — vP' . 
For these models there exists a constant of the motion 

00 

£ = {-n''u^)ng^'^ = Y^'^ p^''< e^^" . (2.65) 
Solving for p leads to 

p = £7g-37/9V-'>'/2 . (2.66) 
This relation is used to eliminate p in the fluid potential 
C/(fl„id) = 24.crx-ie3(2-^)^V-^/2[^-^ _ _ _ 

(2.67) 

There is an additional constant of the motion 113 defined 
by f3 = puz where p = ^n'~^ = -yp^t-^)h ^,|l^. This 
constant of the motion can be used to express y as a 
function of and (3^ . Expressing u^Uq. = — 1 in terms 
of Y and U3 yields the implicit relation 

FrT-^-r + l = 0, (2.68) 



where 

F = ^-2^-2(7-l)„32g2[(37-4),3«+2,3+] (2.69) 

For dust (7 = 1, a case discussed in [|8|) and stiff (7 = 2) 
perfect fiuids this equation is linear and can be explicitly 
solved. For some other values of 7 (namely |, |, |, |, |), 
it reduces to at most a fourth degree polynomial equation 
which can be explicitly solved in principle. However, for 
certain purposes an explicit expression is not required, 
as will become apparent later. Note that the stiff perfect 
fluid potential is the sum of two exponentials. 

The variable 9^ is cyclic and so has a constant conju- 
gate momentum P3 and the equation of motion 

= 122;"ie2"'(/i_)-2p3 . (2.70) 

The single nontrivial supermomcntum constraint requires 

i 

e"' P3 = -2k£v3 , (2.71) 

leading to 

H = ^xriABp^P'' + C/(c) + U^G) + [/(fluid) , (2.72) 
where 

C/(c) = x-'^24:fK'^{v3)^{h^)-^ . (2.73) 

Note that in the type II case and in the type VIq Taublike 
symmetric case (a special case defined by (3~ =0), the 
"centrifugal potential" C/(c) is just an exponential and a 
constant respectively. 

E. Some remarks 

The HH symmetry types include a wide range of mod- 
els describing quite different physical situations. For ex- 
ample, in the diagonal case one may consider any com- 
bination of sources by including the corresponding po- 
tential terms in the Hamiltonian, which together with 
the two possible signs of e leads to numerous space- 
time models that have been considered in the literature. 
Many people have attacked such problems individually 
as though they were completely unrelated to the others, 
each time writing out the field equations and attempt- 
ing to solve them. However, the Hamiltonian approach 
reveals the close mathematical relationship which exists 
between them. 

For example, all the models considered are character- 
ized by a Hamiltonian which can be reduced to the form 

H = ixr;^s/3^/3^ + a;-iC/(taub) - , (2.74) 

where the Taub potential, C/(taub)j is thus just the the 
value of the total potential in the Taub slicing gauge 
X — I. The Lorentz character of the kinetic part of the 
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Hamiltonian suggests that Lorentz transformations 
of the dependent variables play an important role in 
analyzing the dynamics [^o| . Indeed such transforma- 
tions can be used to reveal how many different problems 
are mathematically equivalent. 

In many cases the Taub potential is a sum of exponen- 
tials 



(taub) 



(2.75) 



A Hamiltonian of this type will be referred to as a SE- 
Hamiltonian (for "Sum of Exponentials"). Problems of 
this type are important because the gravitational poten- 
tial and many source potentials are sums of exponentials 
when expressed in the Taub slicing gauge. 

By inspecting the Hamiltonian one can see that some 
of these various problems are either equivalent or closely 
connected. For example, changing the sign of one of the 
variables gaa in the class A gravitational potential ( 2.13 ) 
to go from the SH case to a corresponding static case 
either leaves the Hamiltonian unchanged or is equivalent 
to a change of Bianchi type, e.g., Bianchi types VIII and 
IX are interchanged by this operation. Thus there is an 
isomorphism between various static mod els and the SH 
ones. Furthermore, the MT cases (2.21) have the same 
gravitational potential modulo the sign of ea. When it 
comes to sources the cosmological constant and electro- 
magnetic terms only change sign with the change in sign 
of e. However, a perfect fluid potential differentiates be- 
tween the static and homogeneous cases in an essential 
way since the 4-velocity must be along the timelike di- 
rection in each case. 

The close relationship among many Hamiltonians ex- 
plains the similarity of the expressions resulting from 
solving the field equations for the different models. 



III. THE GENERALIZED FRIEDMANN 
EQUATION 



The generalized Fricdmann equation 
of the form 



a 



i=l 



n| is an equation 



(3.1) 



in a single dependent variable a, where qi are a set of 
distinct constants ordered by increasing value. This gen- 
eralizes the well known Friedmann equation which has 
this form for the scale factor R — e'^. Here the arbitrary 
function Xa allows different choices of the independent 
variable. For convenience, it will be called the slicing 
gauge function. 

The generalized Friedmann equation may be converted 
from exponential potentials to power law potentials by 
introducing the power variable 

w = e'^", (3.2) 



where ^ 7^ is a constant parameter. This leads to 

^2 = x-252^a,u(2+«'/^). (3.3) 

i=l 

A. The power law slicing gauge approach 

To solve the generalized Friedmann equation, it is of- 
ten convenient to introduce a power law slicing gauge 
function lESl 



(3.4) 



where A is a constant parameter. This yields 



w^^^^^a.w'^', r, = 2+|-|. (3.5) 

i—l 

It is convenient to refer to the right hand sides of either 
equation as the "potential" for that variable. 

When only one potential term is present, one can either 
choose A = (71 or A = 2^ = (71 so that a or u respectively 
is affinely related to the independent variable. 

In the case of more than one potential term one can 
always treat one of the terms in the same way one deals 
with the single term in the case with only one potential 
term. However, since there are two parameters available 
(A, (5), one can vary these so that any two power expo- 
nents {ri,rj) assume any pair of real values. Conversely, 
given the values of a pair of the original power exponents 
{qi, qj), the parameters (A, 5) are determined by the val- 
ues of the corresponding new power exponents (r^ , ) in 
the following manner. The scale parameter 5 determines 
the ratio of the power exponent increments 



(3.6) 



after which the additive parameter A is determined by 
either of the relations 



A 



5{n - 2) = q, - 5{r, - 2) 



(3.7) 



In the case of more than one potential term the ob- 
vious way to fix the parameters (A, 5) is to obtain the 
lowest polynomial power variable potential that exists, if 
any. However, even transforming the potential to a poly- 
nomial may be impossible since this requires that the 
original power exponents {qi) be affinely related to a set 
of nonnegative integers. When this is the case for a given 
set of integers, there are two choices of the pair (A, 5) for 
which a polynomial potential occurs, corresponding to a 
positive and negative value of the nonzero power variable 
parameter 5. Thus one always has two different slicing 
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gauges with the same degree polynomial potential, but 
whose dependent variables are related to each other in an 
inverse power relationship. One is an increasing function 
and the other a decreasing function of i? = e" . 

In general, for a case reducible to a polynomial poten- 
tial, if g > denotes the minimum increment between 
the ordered coefficients (qi, . . . , (?„), then the choice 



iA,6)^{qi+2q,q) 



(3.8) 



leads to the polynomial potential with qi corresponding 
to the constant term, while the choice 



(A, S) = (q„ - 2q, ~q) 



(3.9) 



leads to a polynomial potential where the last term as- 
sociated with the final value qn is the constant term. 
These are the minimal degree polynomial potentials that 
are possible. The solutions of the generalized Friedmann 
problem for polynomial potentials up to degree two (four) 
can be expressed in terms of elementary (elliptic) func- 
tions. 

In the case of two potential terms, one can reduce the 
potential to polynomials of first or second degree. The 
parameters (6, A) determining the power variables and 
choices of gauge function which accomplish this are given 
in Table 2, together with the resulting powers. 

If in the case of three potential terms the coefficients 
(91,92,93) are "equally spaced", i.e., 92 - 9i = 93 - 92, 
then there are two values of the pair (A, 5) which map 
them onto either (r^) equal to (0, 1, 2) or (2, 1, 0). Com- 
pleting the square leads to a quadratic two potential term 
problem for a new variable u affinely related to u and the 
same kinds of elementary functions for u result that oc- 
cur for u in two term quadratic case with no first power 
term. For details regarding all of these cases and those 
involving more than three potential terms, see pT|] . 

Often the generalized Friedmann equation occurs as a 
first integral of a second order equation of the general 
form 



for the power variable u = e** when 6^0. Note that 
these have one more exponential potential than the num- 
ber of potential terms in the original second order equa- 
tion. As long as /(e^) has no terms like e~'^^ , then h{e^) 
will also consist of only exponential terms and H(u) of 
powers of u. 

Often the generalized Friedmann equation occurs as a 
decoupled equation in a larger problem. In this case the 
slicing gauge function may be fixed by considerations re- 
lated to the larger problem. When the other part of the 
problem consists of evaluating quadratures, one tries to 
find a slicing gauge so that both the generalized Fried- 
mann equation and the additional quadratures become 
as simple as possible. 



B. The intrinsic slicing gauge approach 

An intrinsic slicing gauge is defined as a slicing gauge 
which relates a linear combination of the variables — 
or an exponential of such a linear combination (a so- 



called "power variable" discussed in section 
affinely to the independent variable, in effect making that 
dependent variable the new independent variable for the 
problem. This type of slicing gauge has played a promi- 
nent role in the context of static HH models, e.g., for 
the spherically symmetric models one usually chooses the 
power variable r as the independent variable. In contrast 
the power law slicing gauges have been much more impor- 
tant for the SH models. These natural selection effects 
seem a bit strange given that the static HH models and 
SH models are so closely related mathematically. 

Recall that the generalized Friedmann equation 
can either be expressed in its "exponential" form 



a = X 

„-2^2 



or its 



x^'O'^^^^^QiU^^'^'^ where u = e 
intrinsic slicing gauge means that one either chooses a 
or some u variable as the independent variable. Such 
choices are accomplished by setting 



"power" form u — 
In this context an 



+ S9' + fie'')=0, 



(3.10) 



where f{9) is a sum of exponential terms. This equation 
has the first integral 



2S0 rn 2 



h{e')] , hie") = -2e-''' / e'''f{e') dO , 

(3.11) 



which may be rewritten in the form of the generalized 
Friedmann equation as 



or in the power form 



(3.12) 



u 



H{u) + S^£, H{u)^5'^u^h{u^'^) (3.13) 



n 

= A: x„ = (5(^a,u2+9-^)i/2. 



(3.14) 



These choices may seem trivial. However, if the gen- 
eralized Friedmann equation is part of a larger problem, 
they may not be. For example one often encounters prob- 
lems where one has a cyclic variable which is given by 
a quadrature of the form (3^ — J x~^ri^^pBdX, where 
Pa is the constant conjugate momenta associated with 
(3^. Choosing a power variable u leads to the integral 

pA ^ rj'^Bp^g J(^^^^^^^2+g,5)-l/2 ^j^j^j^ gj^gg ^j^g 

dependent variable as a function of the independent one. 
If on the other hand one chooses a power law slicing 
gauge, this may simplify the equation for the /3"^ vari- 
able while transfering the difficulties to the generalized 
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Friedmann equation. Even if this latter equation is solv- 
able in terms of a quadrature, it gives the independent 
variable as a function of the dependent variable u, in 
contrast with the intrinsic approach which yields the de- 
pendent variable as a function of the independent one. 
Thus one is faced with the problem of trying to invert a 
quadrature in the power law gauge approach, something 
which often fails to lead to familiar functions. Thus the 
intrinsic approach sometimes has advantages. As an ex- 
ample, see p^ ], where the intrinsic approach was used 
in the context of SH Bianchi type V orthogonal perfect 
fluid models. 



IV. KILLING TENSOR SYMMETRIES AND 
HOW TO USE THEM 

To solve a Hamiltonian problem, one needs to find and 
exploit symmetries which lead to constants of the mo- 
tion. In the search for such symmetries, a particular 
slicing gauge will be introduced which makes the equa- 
tions of motion equivalent to the geodesic equations as- 
sociated with a certain metric (not to be confused with 
the spacetime metric). This makes it natural to look 
for a particular type of symmetry, called a Killing tensor 
symmetry, which generalizes symmetries corresponding 
to cyclic variables and Hamilton- Jacobi separability. 

However, this slicing gauge is not usually well suited 
to exploiting the symmetry so that explicit exact solu- 
tions can be obtained. The existence of constants of the 
motion is not sufficient to obtain such solutions explic- 
itly. This requires the stronger condition of a decoupling 
of the equations of motion (which in turn leads to con- 
stants of the motion), for which other slicing gauges are 
needed. Even within this latter class of slicing gauges 
there are choices of gauge which lead to simpler forms of 
the exact solutions. The situation is similar to the case 
of a cyclic variable. Any slicing gauge which doesn't in- 
volve this cyclic variable leads to a constant of the motion 
and decoupling. However, some choices lead to simpler 
equations and simpler expressions for their solutions. 

The present section discusses the symmetry properties 
of Hamiltonians and how to explicitly find solutions of the 
equations of motion in as simple a form as possible. It 
then relates these results to the Hamiltonians of the type 
encountered when dealing with the HH models described 
in section 2. 



A. Why Killing tensor symmetries? 

In classical mechanics one usually encounters Hamilto- 
nians of the form 



H 



U ^E, 



(4.1) 



where the symmetric matrix (g°^) is positive-definite. 
For such problems there exists an elegant geometric re- 
formulation of the corresponding equations to a geodesic 



flow on a certain geometry ||2^,|26|,^ . To accomplish such 
a reformulation one first introduces a new Hamiltonian 
TC = H — E = 0. A Hamiltonian system of this kind can 
be reparametrized by choosing a new independent vari- 
able A in place of the original one A, leading to a new 
Hamiltonian 



Hm = Mn= N{H - E) , Af=d\/d\. 



(4.2) 



The final step is to make the Hamiltonian purely ki- 
netic. This is accomplished by a particular choice of 
parametrization, Mj = ^{E — U)^^ , and by adding a 
constant to the Hamiltonian: 



Hi 



1 



4(S - U) 



g^'paPb 



(4.3) 



The corresponding Lagrangian equations can then be 
reinterpreted as those of the geodesic fiow of the so-called 
Jacobi metric Jab = 2{E — U)gab, where gab is the matrix 
inverse of 5°^, and A is an affine parameter along each 
geodesic. 

The geodesic reformulation does not rely on the 
positive-definiteness of gab- It works locally for any non- 
degenerate indefinite matrix gab- For the HH spacetimes 
of section 2, the Hamiltonian can be put in the form 



H 



(taub) 



0. 



(4.4) 



The Hamiltonian is already parametrized and must van- 
ish so there is no need for the first steps in the above 
procedure. The choice 



X^Xj ^ 2 1 [/(taub) I 



(4.5) 



leads to a Hamiltonian with a constant value of the poten- 
tial. Redefining the Hamiltonian by this constant yields 



Hj = Tj = \JabI3^P'' = -isgn(C/( 



taub) } 



(4.6) 



where Jab = xj rjAB - In contrast with the usual classical 
mechanical problems, the underlying geometry here is 
Lorentzian rather than Riemannian. 

To solve a Hamiltonian problem one needs to find 
enough symmetries leading to constants of the motion, 
i.e., variational symmetries, which can be used to reduce 
the problem sufficiently so that the reduced problem can 
be solved. The Jacobi formulation is especially suitable 
for finding such symmetries since they and their associ- 
ated constants of the motion take a particularly simple 
form in this formulation as discussed below. Further- 
more, the geometric framework makes available a wide 
range of tools familiar from symmetry investigations on 
ordinary spacetime. 

Variational symmetries are transformations of the 
phase space which can be represented as transformations 
on the tangent bundle (velocity phase space) generated 
by vector fields of the following form |28| 



d 

dx"- 



(4.7) 
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The simplest variational symmetries are the point sym- 
metries of the configuration space itself. In the Jacobi 
formulation, a generator (j)"" {x)d / dx"" of a point symme- 
try corresponds to a Killing vector field of the Jacobi met- 
ric ]29| ]. All other variational symmetries involve deriva- 
tives of the dependent variables and are called generalized 
symmetries |2^. 

The simplest generalized symmetries are the ones for 
which the components of their generating vector fields 
are linear and homogeneous in the derivatives, i.e., (j)"" = 
K°-b{x)x^ . This corresponds exactly to a second rank 
Killing tensor K"'}, of the Jacobi metric |2^. Recall that 
a Killing tensor Kab is a symmetric tensor for which 
K{ab;c) — 0; this includcs the trivial case in which Kab 
is proportional to Jab and the special case in which 
the Killing tensor is the symmetrized tensor product of 
Killing vectors. Killing vector symmetries give rise to 
constants of the motion which are linear and homoge- 
neous in the momenta, while second rank Killing tensor 
symmetries correspond to homogeneous quadratic con- 
stants of the motion {K°-^paPb = const.) |^,|oj. This is 
another way of understanding how Killing tensor sym- 
metries are the natural generalizations of Killing vector 
symmetries. 

The separability condition for the Hamilton- Jacobi 
equation is equivalent to the existence of a sufficient num- 
ber of second rank Killing tensor symmetries of the Ja- 
cobi metric Jab pOj . Furthermore, in the case of an indef- 
inite Jacobi metric, separability requires the Killing ten- 
sors to have nonnuU cigendirections. However, it turns 
out that Killing tensors with null eigendirections can 
also be used to solve the Hamiltonian equations. Thus, 
in contrast with the positive-definite case. Killing ten- 
sor symmetries in the indefinite case are more powerful 
than Hamilton- Jacobi separability when searching for ex- 
act solutions. 

The Jacobi slicing gauge is not the most convenient 
gauge to use when it comes to actually producing ex- 
act solutions. There are more suitable symmetry com- 
patible choices of independent variable. However, other 
such choices lead to corresponding constants of the mo- 
tion which are quadratic but not homogeneous in the 
momenta. Since the constant of the motion takes an es- 
pecially simple form in the Jacobi slicing gauge, entirely 
characterized by the Killing tensor alone, it is this gauge 
which makes finding the associated Killing tensor sym- 
metries as simple as possible. Thus it seems clear that 
a search for these symmetries is a natural first step in 
attempting to find exactly solvable problems. 

Many of the HH field equations have long been known 
to admit a Hamiltonian formulation very much like the 
traditional classical mechanical problem except for the 
signature difference. Therefore it is surprising that no 
attempt has ever been made to solve the field equations 
by means of Hamilton-Jacobi separation given its success 
in the Riemannian case. It is perhaps more understand- 
able that no one has used the more powerful but less 
familiar Killing tensor techniques. 



B. Killing tensor symmetries 

For many of the known exact HH solutions, the field 
equations are ultimately expressible in terms of two non- 
trivial degrees of freedom (although more variables may 
be involved). It is therefore useful to start with the case 
of two de gree s of freedom, reviewing and extending ear- 
lier work [ pi[ . 

When dealing with 2-dimensional geometries the group 
of conformal transformations plays a particularly impor- 
tant role since in this case it is infinite-dimensional. For 
2-dimensional Lorentzian geometries it is useful to use 
null variables since they are closely related to this group. 
A Lorentzian 2-metric can always be written in the form 

ds^ ^ -2G{w,v)dwdv (4.8) 

in terms of the null variables w and v. A general con- 
formal transformation to new null variables V and W is 
of the form v — F{y),w — F{W), and results in a new 
conformal factor G = F'{V)F'{W)G. One may then 
introduce new conformally inertial coordinates which di- 
agonalize the metric 

T ^\{W + V), X^l{W~V). (4.9) 

It turns out to be convenient to classify the geometries 
allowing Killing tensors into three cases characterized by 
different conditions on the form to which the conformal 
factor can be transformed. These three cases can be char- 
acterized as follows pH : 

(i) (null) : G = [A{W)V{V) + B{W)][dV{V)ldV] , 
(ii) (nonnuU H-J) : G = G{T) + D{X) , 
(Hi) (nonnuU harmonic) : 

G yy -I- G = i(G,TT + G^xx) = . (4.10) 

In the above expressions A, B, G, D, V are arbitrary func- 
tions and the variables V, W or T,X will be referred to 
as symmetry-adapted variables. 

The Killing tensor in case (?) is characterized by hav- 
ing a degenerate eigenvalue corresponding to a single 
null eigenvector. The remaining cases {ii) and {Hi) have 
a Killing tensor characterized by nonnuU eigenvectors. 
The case {Hi) corresponds to a conformal factor satis- 
fying Laplace's equation which is therefore a harmonic 
function of W,V or equivalently of T,X. Of the three 
Killing tensor cases it is only the case {ii) that correponds 
to Hamilton-Jacobi separability. Cases (i) and {ii) are of 
special importance. The case («) will be referred to as 
the "null case" while, for simplicity, the case {ii) will be 
referred to as the "nonnuU case" unless explicitly quali- 
fied as the "H-J case" . The case {Hi) will be referred to 
as the "harmonic case." 

The expression for the conformal factor in case {i) can 
be simplified to 
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(ib) G ^ A{W)V + B{W) 



(4.11) 



by changing the variables to V — J V,y dV and W — W. 
However, the original expression is more convenient to 
use as starting point for an analysis of some of the ex- 
amples to be discussed below. Note that the form {i) or 
{ih) is invariant under the transformations — > VF(H^), 
leading to an equivalence class of symmetry- adapted de- 
pendent variables. 

When a Killing v ector exists rather than a nontrivial 
Killing tensor, Eq. ( 4.1Cl| ) continues to hold with either 
AB = or CD = 0. The condition AB = is associated 
with the existence of a null Killing vector, leading to a 
flat Jacobi geometry. The condition CD = Q corresponds 
to a nonnuU Killing vector field, with C = describing 
the timelike case and D = the spacelike case. 

The case (u) is also generalizable directly to a 2- 
dimensional subblock of the metric in higher dimension 
or to the case of complete separability where this factor 
has the form 



G = C(T)+^A(^^). 



(4.12) 



The eigenvectors of the Killing tensor in the subblock 
case are aligned with the conformally inertial coordinate 
directions T and X. In the case of complete separability 
one has a set of nonnuU Killing tensors whose eigenvec- 
tors are aligned with T and X*. Intermediate cases be- 
tween these two extreme cases are also possible. In the 
special case of a Killing vector, one can choose symmetry- 
adapted conformally inertial coordinates such that the 
conformal factor is independent of one of them. 



C. How to use Killing tensor symmetries 

We are not just interested in obtaining a sufficient 
number of constants of the motion to solve the problem 
implicitly. We want to obtain explicit exact solutions. To 
do this one needs to decouple a sufficient number of equa- 
tions. In general the Jacobi slicing gauge does not allow 
this even if one has a sufficient number of Killing ten- 
sor symmetries. However, there are other choices of the 
independent variable which permits this if one chooses 
symmetry- adapted dependent variables. Since there is 
some freedom in the choice of the dependent variables, 
this freedom may be exploited to help simplify the prob- 
lem further. A specific choice of dependent variables of- 
ten suggests which independent variable to use to actu- 
ally integrate the equations of motion. 

To simplify the discussion only 2-dimensional examples 
will be considered. To find useful slicing gauges one must 
first re-express the general-slicing-gauge Hamiltonian in 
terms of the symmetry-adapted dependent variables and 
the special conformal factor G associated with the ex- 
istence of one of the three types of Killing symmetries. 
This then suggests a choice of independent variable in 



terms of which the equations of motion can be explicitly 
solved. 

Suppose one has a Hamiltonian of the form 



H = \x{-dP + /?2) + a;-iC/(taub) ■ 
By introducing the null variables 

w = a + P ^ V = a — f3 , 
one transforms the Hamiltonian to the form 



H 



(taub) 



(4.13) 
(4.14) 
(4.15) 



This Hamiltonian corresponds to a Jacobi geometry with 
the conformal factor G ~ |C/(taub)| (recall that xj = 

2|C^(taub)l)- 

Re-expressing the Hamiltonian in terms of a new set 
of null variables w — wiW) and v — v{V) leads to 



H 



^x{dw/dW)idv/dV)WV + x-'^U^t^^,bj . (4.16) 



If the Jacobi geometry admits a Killing tensor and these 
new variables are associated symmetry-adapted vari- 
ables, then one can make the identification 

G - idw/dW)idv/dV) I [/(taub) I . (4.17) 

Inserting this relationship into the Hamiltonian leads to 

H=^-^x I [/(taub)!"' GWV + a;-iC/(taub) • (4.18) 



It is convenient to introduce a new slicing gauge func- 
tion y defined by 



y = |C^(taub)l ' G: 



(4.19) 



when discussing symmetry-adapted slicing gauges. Ex- 
pressing the above Hamiltonian in terms of G and y yields 



H 



yirF + y-isgn(C/(taub))G. 



(4.20) 



The two cases (z) and (m) in Eq. ( [4.10 ) will be treated 
separately when it comes to finding slicing gauges which 
lead to decoupling and explicit solution of the equations 
of motion. The case (Hi) does not seem to have the same 
importance as the first two Killing tensor cases and will 
only be commented on briefly when it does occur. 



1. The null Killing tensor case 

The Lagrangian equatio ns of motion corresponding to 
the previous Hamiltonian (4.20) are 



W + y'\dy/dW)W^ 

-2y-'d{y-' sgn(C/(taub))G)/9V - , 
V + y-'^{dy/dV)V^ 

-2y-'d{y-' sgn(C/(taub))G)/9Ty = . (4.21) 
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Assume that one has chosen symmetry-adapted variables 
so that the conformal factor takes the simpler form (ib) 



sgn(f/(taub))G = A{W)V + B{W) . 



(4.22) 



It is then possible to choose a slicing gauge so that the 
equation of motion for the variable W decouples. Those 
slicing gauges which allow this decoupling are character- 
ized by the condition 



y = y{w) 



(4.23) 



and will be referred to as decoupling slicing gauges. A 
decoupling of the variable V occurs for the analogous 
conditions with V and W interchanged. 

For a given choice of y{W) one can choose new depen- 
dent variables W and V defined by 



dW/dW = y{W) . 



V ^V. 



(4.24) 



This transformation eliminates the quadratic first deriva- 
tive term from the decoupled equation and leads to the 
Hamiltonian 

H = -\WV + \dW{W)ldW\\A{W{W))V + B{yV{W))\ 
= -\WY + A{W)V + B{W) , (4.25) 

which has its kinetic energy in the standard Minkowski 
null form. 

A Hamiltonian of this final form leads to the decoupled 
equation 



W = 2A{W) , 
which has the first integral 



(4.26) 



(4.27) 



= AJ A{W) dW + const . 
Rewriting this in the form 

^W^+U{W)^E, U{W) ^ -2 [ A{W)dW , (4.28) 



where E is a. constant, allows U{W) to be interpreted as 
a potential for this 1-dimensional problem. 

By appropriately choosing the function y(W), one can 
make the function A(iy) take any desired form. How- 
ever, there are preferred choices. For example, choosing 
A(VF) to be proportional to one of the powers 0, 1, or 
—3 of or to a sum of terms involving the powers and 
1 of W leads to a generalized Friedmann problem with 
elementary function solutions for W as a. function of the 
independent variable. The first two single term choices 
lead to a linear and quadratic potential, while in the last 
single term case the transformation W — leads to 
a linear potential. The combination of and 1 power 
terms leads to a quadratic polynomial potential. All of 
these potentials were discussed in section [H. After a 
choice of gauge, leading to a particular form for ^(M^), 



the remaining variable V is determined by first solving 
the Hamiltonian constraint for V, yielding 

V ^ (2/W)[A(W)V + B{W)] = WV/W + 2B{W)/W , 

(4.29) 

and then integrating this linear first order equation, lead- 
ing to 



V = W [ const + 2 / B{Wy^dX 



(4.30) 



Whether or not the integral can be expressed in terms of 
elementary functions depends on the specific case one is 
dealing with. 

It follows from the above formulas that the expressions 
for C/(taub) X then become 

f/(taub) = [dW{w)/dw] [dViv)/dv] [AiW)V + B{W)] 
= [g{w)V{v) + h{w)][dV{v)/dv] 
x^[dW{w)/dw][dV{v)/dv]. (4.31) 



As noted above after Eq. (4.11), the Jacobi geometry is 
fiat if g{w) = or h{w) = 0. One can often see by 
inspection whether or not a Taub potential is of the above 
form. For such cases it is convenient to start with the 



gauge function x and derive Eq. ( [4.25]) . Those shcing 
gauges which allow decoupling are characterized by the 
condition 



X — Xi{w) dV{v)/dv , 



(4.32) 



where xi{w) takes the role y previously played. Decou- 
pling in the Taub slicing gauge itself can only occur if 
V is a linear function of w. If 1/^ is not a linear function 
of V one makes a variable change from w to 1^ to obtain 
decoupling. This leads to 



H ^ -^Xi{w)wV + xi {w)-^[g{w)V + h{w)] ^ , 

For a given choice of xi{w), the choice 

dW/dw — xi{w) , 

of a new dependent variable W leads to 

X = [dW{w)/dw] [dV{v)/dv] , 
H = -\WV + A{W)V + B{W) , 



(4.33) 



(4.34) 



(4.35) 



that is, the same expressions as the end result of taking 
the "y approach" . 

Many Hamiltonians correspond to null Killing tensor 
cases, thereby explaining the existence of a whole host of 
exact solutions. 
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2. The nonnull Killing tensor case 



Making the transformation T = ^{W + V), X 



\{W — V) transforms the Hamiltonian in Eq. ( 4.20| ) to 



H^\y {-f^ + X^) + sgn(;7(taub)) G . (4.36) 

If one has a nonnuU KiUing tensor case expressed in 
symmetry-adapted variables so that 



X = |C/(taub)|G ^ 



sgn(C/(taub))G' = C(T) + i^(X) 
then the choice 

y=l 

leads to the Hamiltonian 

H = i(-r2 + X^) + C{T) + D{X) , 
for which both the T and X equations decouple 

f ^ dC{T)/dT , X ^ -dD{X)/dX . 
These equations give the first integrals 

-if^ + CiT)^Et, ^X^ + D{X) = E,, 
where the integration constants have to satisfy 



(4.37) 
(4.38) 
(4.39) 
(4.40) 
(4.41) 
(4.42) 



because of the Hamiltonian constraint. Equation ( 4.39| ) 
corresponds to Hamilton-Jacobi separability and is gen- 
eralizable to any dimension. 



D. Power variables and power law slicing gauges 

Power variables are often the simplest choice of de- 
pendent variables adapted to Killing tensor symmetries. 
Such a variable is defined by an expression of the form 



HQ) 



— "(0)6 — "(0)6 



(4.43) 



Similarly power law slicing gauges are often the simplest 
slicing gauges which lead to decoupling. Such gauges are 
characterized by the following form of the slicing gauge 
functions [Esl 



(4.44) 



The conditions Qi — Q2 — Q3, or equivalently Q± = 
or Qx = as appropriate, (and (54 = if a scalar field 
is present) characterize the isotropic power law gauges 
which can be useful for SH models in the case of spatially 
isotropic sources like an orthogonal perfect fluid or a cos- 
mological constant term, the Taub time gauge being an 



obvious example. However, because of anisotropic spatial 
curvatures and sources, one needs to exploit anisotropic 
dependence of the function N on the individual scale fac- 
tors to make progress in obtaining solutions or simplify- 
ing the equations. The condition Qi + Q2 + Q3 = 1 or 
equivalently Qo = 1 leads to a scale invariant indepen- 
dent variable. The conformal gauge corresponds to both 
conditions holding, leading to N = Af(o)e^ . Finally it is 
worth noticing that power variables and power law slic- 
ing gauges are not only useful when it comes to finding 
symmetry-adapted variables but also in the context of 
qualitative analysis. 



E. Killing tensor symmetries in models 
characterized by SE-Hamiltonians 

1. Null decoupling 

Consider a potential of the form 

c/(taub) = [giy^Wiv) + hiw)][dv{v)/dv] . 



(4.45) 



If g{w), h{w) and V are sums of exponentials including 
possible constant additive terms, or if g{w) — and V is 
linear in v, then one has a SE-Hamiltonian corresponding 
to a null Killing tensor case. Of particular interest is the 
case in which is a single exponential term or linear 
in V in the "flat" g{'w) = case. The remainder of the 
subsection analyzes this situation. The Taub potential 
for the cases 

y = e^«7co, co^O, giw)^0 or giw) = , 
V = v, co = 0, giw) = 0, (4.46) 

can then be formally written as 
C/(taub) = [coiJ2 Are''ne'°" + E B^e'^""V'''" ■ (4.47) 



i=0 



i=0 



Note that the three separate cases cq = 0, or = or 
Bi = for all values of i all correspond to fiat cases. 

Choosing x = xi{w)e'^°^ and V = e'^"^/co as a new 
variable in the cq 7^ case leads to the Hamiltonian 



H 



(4.48) 



In the case cq = 0, C/(taub) is independent of v, and the 
Hamiltonian takes this same form with V equal to v. 

Usually the simplest class of choices for the remaining 
gauge function xi{w) are power law slicing gauges. For 
such gauges this function is proportional to an exponen- 
tial, leading to 



X = Be' 



(4.49) 
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and the Hamihonian 
i7 = -iBe"™wF 



(4.50) 



i=0 1=0 

Here one may choose m to be either zero or nonzero. 

a. The gauge choice m — 0. Setting B = 1 one then 
has 



H = -^wV+ [co2( Ae'^nV + j2 ^^e' 



(4.51) 



This Hamiltonian leads to a decoupled second order dif- 
ferential equation which in turn can be integrated to yield 
a first order differential equation 



4co^y(A,e"'™)/a, + const 



(4.52) 



c. The case Ai — for all i (the flat case). 

The gauge choice m = 0. This leads to w = 0, 
so that w is affinely related to the independent variable. 
Solving the Hamiltonian constraint for V leads to an eas- 
ily integrated expression for V (or w if cq = 0) which re- 
sults in a linear combination of exponentials unless one 
of the exponential coefficients bi vanishes, leading to a 
linear term. 

The gauge choice rn ^ 0. This leads to = and 
W is affinely related to A. Setting the additive constant 
to zero, again one can integrate the expression for V to 
obtain T^, this time resulting in a linear combination of 
powers of the independent variable 



hi jr. 



(4.55) 



if all the coefficients are nonzero. If not, the single 
Oi — Q term gives rise to a term linear in w which has to 
be added to this expression in order to give the correct 
form for the first order equation. 

The decoupled first order equation for w arising from 
the Oi ^ case is simply the generalized Friedmann 
equation, while V may be obtained formally by solving 
the Hamiltonian constraint (which is linear in y) for its 
derivative and integrating. To simplify the process of in- 
tegrating this final pair of equations (abandoning further 
use of the Hamiltonian/Lagrangian approach), one can 
reintroduce the gauge freedom that was fixed above by 
a specific choice of x\(yj) made to obtain the decoupled 
first order equation. However, usually this extra freedom 
is not needed to find many simple explicit expressions for 
the solution of the system. 

6. The gauge choice m ^ 0. If one instead chooses 
m ^ 0, and follows the procedure after Eq. ( 4.33| ) and 
uses the power variable 

W = (4.53) 



and sets B ~ m, then one obtains 
H = -\WV 



(4.54) 



6i/m — li 



In the case that all are nonzero, the decoupled equation 
gives rise to the power form of the generalized Friedmann 
equation. An a.; = term yields a logarithmic term in 
the integrated expression of the decoupled second order 
differential equation. As in the previous case, to integrate 
this final pair of equations, one can reintroduce the slic- 
ing gauge freedom involving rescaling of A derivatives by 
functions of W in order to simplify the problem. With- 
out resorting to this additional freedom, one can make 
the following case by case analysis of some types of SE- 
Hamiltonians which often arise. 



unless one of the coefficients hi vanishes, leading to a log- 
arithmic term. (Note that if cq = all the above expres- 
sions for the Ai = Q case remain valid with V replaced 
with V.) 

d. The case co ^ 0, ^ 0, A, = (i>l ). 

Subcase oq 0. Here the choice to = is appro- 
priate, leading to iu ~ const and quadratic solutions for 



Subcase ao ^ 0. Here the three choices m = 
ao,±^ao lead to element ary fu nction solutions for W as 
described after equation (4.28). 



e. Two examples. 

One null exponential potential term in any di- 
mension. When a problem is characterized by a single 
null potential term, one first uses a Lorentz transforma- 
tion from the original conformally inertial coordinates to 
a new set l3°,f3,(3^ so that the Hamiltonian takes the 
form 

H = + + ^'^] + \x-'Ke^'^'^"+P^ . 

s 

(4.56) 

Then one introduces the null variables 

w = (3° + (3^ v = p^ -p. (4.57) 
In these variables the Hamiltonian then takes the form 



H = —-^xv w 



^x 



[0^ 



Ke 



(4.58) 



where — "^gps^ arises from nonnuU cyclic variables 
(3^ and S denotes the index labeling the nonnuU cyclic 
variables, assuming that x is chosen to be independent 
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of them. The variables (5^ can be evaluated in terms of 
quadratures arising from the equations 



(4.59) 



The full problem can therefore be decomposed into a 
set of quadratures and a 2-dimensional problem coming 
from the above reduced Hamiltonian for w and v. This 
is a flat null Killing tensor case, Ai = 0, which has been 
discussed above. One need only identify the constants 
Co = 0, 6o = 0, bi = C, Bq = -^O^, Bi = jK in equation 



(4.48) 



One exponential potential term in two dimen- 
sions. The one potential term case corresponds to a 
flat Jacobi geometry. The single potential term can be 
identified in two ways with the explicit terms in equation 
( [4.47D . Either AqCq ^ or Bq =^ 0. These two different 
identifications lead to different natural slicing gauges. 

If Aq — 0, then Bo ^ 0, then one can choose m — and 
a translation in A leading to w X. Integration of the 
Hamiltonian constraint yields an exponential expression 
for V unless 6o = in which case V ~ A + const. In 
this latter case the Taub slicing gauge potential is null, 
i.e., it depends only on the null variable v. If one chooses 
instead m ^ 0, then one can make a translation in A 
such that - A and y - A''"/™ + const if bo ^ 
or F ^ InA + const if bo = 0. In the first case the 
Taub potential is nonnuU and the choice m = bo makes V 
affinely related to the independent variable. When both 
variables are affinely related to A, the choice of slicing 
gauge makes the single potential term in the Hamiltonian 
a constant, i.e., the gauge is the Jacobi slicing gauge. 

If AoCo ^ 0, then either ao = or oq 7^ 0. If ap = 0, 
then with m = 0, the variable w is quadratic in A and the 
expression for 'In V is a standard integral. If oq ^ then 
the three choices m = ao,±iao are relev ant according 
to the discussion following equation ( 4.27 ). The choice 
m = ao leads to a quadratic solution for W. Choosing 
the zero of the independent variable to eliminate the lin- 
ear coefficient leads to power law solutions for V, i.e., 
V ^ X"' , where a is determined by the quadratic coef- 
ficient and B. For the choice m = ^ao, W is a hyper- 
bolic/trigonometric sine or cosine and the expression for 
'\nV' is a standard integral. For the choice m — — ^oq, 
W'^ is quadratic in A leading to 'InV^ being a quartic 
expression. An entirely parallel discussion holds with w 
and V interchanged corresponding to those cases where v 
decouples. This will lead to additional slicing gauges. 

/. Linear decoupling. It is not necessary to choose 
null variables to solve problems admitting null Killing 
tensors. The only requirement is that a variable and 
its derivative occur linearly in the Hamiltonian. As an 
example, the Hamilton ian p roblem corresponding to the 
potential of equation ( 4.47 ) can be treated as follows. 
Making a transformation 



V = {In g + w)/ CO ^ g ^ 6"°''-'" , (4.60) 
and choosing a new gauge function y ~ x/ g leads to 



H = -y{—){wg + gw'^) 
2co 



(4.61) 



i=0 1=0 

This Hamiltonian leads to the decoupled equation 
w + (2/-1 ^ - l)w^ - 2co'y-' J2 A^e^'"^'^- = . 



(4.62) 



For the choice y = De'^"' this leads to 
2co 



(d _ l)^/;2 - y A,e('^'-2d+2)^„ ^ ^ (4 



which is just the second order form for the generalized 



Friedmann equation discussed in section HI 



Note that this discussion starting from Eq. (4.61) easily 
generalizes to include all real values of g. The identifica- 
tion y = — 4coe^'"z shows that the Hamiltonian (2.45) of 
the Taub-NUT-M spacetime is of this form. 



2. NonnuU decoupling 

A particular type of nonnuU decoupling occurs for SE- 
Hamiltonians in any dimension when after a Lorentz 
transformation to new conformally inertial coordinates 
(3^ (where (3^ denotes the timelike variable), each expo- 
nential potential term involves only a single new (non- 
null) variable. In other words the Taub potential takes 
the form 



(taub) 



iAe 



(4.64) 



Then the total Hamiltonian is just the sum of indepen- 
dent Hamiltonians Ha for each new variable, constrained 
only by the Hamiltonian constraint on the sum of the in- 
dividual energies 



HA^^VABiP'^y + Y.B.Ae'' " ^Ea, (4.65) 



Each equation Ha — Ea is a generalized Friedmann 
equation. 

The simplest example of this occurs for a single (non- 
null) potential term. The Hamiltonian expressed in terms 
of the new variables then has the form 
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H = ^xtjabP-^P^ + ^x^^Ke^''^'' = (no sum over P) 

(4.66) 

In the timehke case the index P will assume the value 
0, while in the spacelike case it will assume one of the 
remaining allowed values. 

Letting x only depend on (3^ results in the equations of 
motion [x^^)' = for the remaining cyclic variables 
Q ^ P. These lead to constant values of the momenta 
pQ — xrjQuf]^ . Thus the cyclic variables are determined 
by the equations 



(4.67) 



The noncyclic variable may be solved for using the Hamil- 
tonian constraint, which may conveniently be rewritten 
in the standard form of a generalized Friedmann problem 
with two exponential terms 



~2[i;2 + Ke'^''^''] , P = ^ timelike case , 
-2[r - Ke^'''^''] , P^O^ spacelike case , 

(4.68) 

where — J2qPQ^ ^^'^ ^ — ~'n^^PQpR- This problem 
is easily solved in various choices of slicing gauge and 



power variables as discussed in section [II 



There are 2-dimensional SE-Hamiltonians correspond- 
ing to nonnull Killing tensor cases which require nontriv- 
ial conformal transform ation s as well. If one starts with 
the Hamiltonian in Eq. ( 4.36| ) and assumes that the con- 
formal factor is a polynomial in T and X of the following 
form 



sgn(t/(taub)) G = ^(A,r + B,X') , 



(4.69) 



i=0 



which corresponds to a simply solved separable problem, 
then the variable transformation 



X = e° 



(4.70) 

where a and b are arbitrary nonzero constants, and the 
relation (following from the definition (4.1£)) 

y = {Aabe'''"+'"'y^x 

leads to the Hamiltonian 



(4.71) 



i=0 /c=0 
^{i-k+l)aw+(k+l)bv 



Thus the identification 



(taub) 



(4.72) 



(4.73) 



(4a6) ^ ^(A, + i^lfB,) f ^ J e[^-k+i)a^+(k+i)bv 



i=0 k=0 



can be made. There are several solvable cases appear- 
ing in the literature corresponding to polynomials of low 
degree in Eq. (4.69). 



3. Nonnull-null decoupling 

For those cases whose Jacobi geometry allows both a 
null Killing tensor and a nonnull Killing tensor, one may 
take yet another approach in solving the field equations. 
For example, if the Taub potential of a 2-dimensional 
problem is both of th e form (4.47), corresponding to w 
decoupling, and ( 4.64 ), then the Hamiltonian must take 
the form 



H = -^XWV + X ^[/(taub) , 



(4.74) 



A similar expession holds for v decoupling with w and v 
interchanged. 

One may then boost (w^v) = {kw,k^^v) to {w^v) — 
(a -I- (3, a — (3), where fc > is determined so that each 
exponential term in the potential depends on only one of 
the new conformally inertial coordinates a or (3, namely 
by the condition 



C = 2&fc = ±2cA:^\ ^ fc=|c/6|i/2 
Then if sgn fee = 1 , the potential becomes 



(taub) 



Dab 



-C/3 



(4.75) 



(4.76) 



while if sgn 6c = — 1, a and /3 are interchanged in the 
potential. 

Assuming sgn 6c = 1, then two power law slicing gauge 
choices lead to mutual decoupling of a pair of variables, 
the null variable w and one of the new nonnull inertial 
coordinates. These correspond to making the D2 and D4 
potential terms constant respectively. These choices are 



case (a) : 
case (6) : 



X 
X 



e^" ; w and a decouple, 
w and P decouple. 



-CP . 



(4.77) 



with the respective Hamiltonians 

case (a) : H = -^e'^''w{2a - w) + Die'^'' + D2 



case (6) 



H 



-ie-^^zi;(7i) - 2/3) + Die'^^^^-'^ + D^e^"" 



(4.78) 



Similar results hold with a and f3 interchanged for 
sgn 6c = — 1 . 

The decoupled equations for case (a) are 



0^ SL/Sa = e^^i-ijj- ^Cw-" 
= 5L/Sw + 5L/5a + CH 



-2CtSi 



(4.79) 



while those for case (6) are 
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{) = 6L/5I3=- 







5L/Sw ■ 



-SL/SP- 



CH 



C/31 



(4.80) 



Each of these four deco upled equations is equivalent 
to the second o rder form ( j3.10| ) of a generalized Fried- 
mann equation ( 3.12| ). They have the following respec- 
tive values of the potential term /(e^) and the constant 
6 appearing in the latter equations 



fie') 



(4.81) 



The equivalent generalized Friedmann equations are then 
respectively 



6|2 = h{e'')+£e 



-258 



Di + D2e-'^" + £e^2Ca 



^2^ 

he' 



(4.82) 



Finally the following power variables convert these equa- 
tions into quadratic potential problems 



D2U + S, 
D3+£u, 
D4U + £ . 



u 
u 



(4.83) 



The solutions of these equations describe 1-dimensional 
motion in a quadratic potential and lead to solutions for 
the dependent variable u which are afhnely related to ex- 
ponential or trigonometric or hyperbolic sines and cosines 
of an argument afhnely related to the independent vari- 
able. Of course this method can also be used to treat the 
single potential term case, but the previous two methods 
are simpler. 



4. Power variables, power law slicing gauges and 
SE-Hamiltomans 

If one looks at the literature on exact solutions one 
almost always find the solution expressed in power vari- 
ables and power law slicing gauges. Why is this the 
case? The answer is that practically all solvable prob- 
lems are described by SE-Hamiltonians with a relatively 
few number of potential terms. As seen above, the sim- 
plest symmetry-adapted variables and slicing gauges are 
usually power variables and power slicing gauges directly 
related to these exponential terms. On the other hand, if 



one has many exponential terms one might be forced to 
use non-power law variables. For example, this happens 
when the function V, in the null Killing tensor case (i), 
is a sum of exponential terms. 



F. The intrinsic approach to null Killing tensor 
problems 

Apart from the general case of the Jacobi slicing gauge, 
the only specific slicing gauges which have been consid- 
ered here are the power law gauges. This subsection will 
show how another important class of slicing gauges arises 
in a natural way for the null Killing tensor cases. These 
gauges are the so-called "intrinsic" slicing gauges. 

Suppose one chooses symmetry-adapted variables and 
a symmetry compatible slicing gauge in the null Killing 
tensor case so that one obtains the Hamiltonian (4.54) 



H = -\WV + A{W)V + B{W) = . 



(4.84) 



Then one obtains the first order decoupled equation for 
W 



W'^ = 2[E^U{W)] = F{W) , 
where U{W) ^ -2 



A{W)dW . (4.85) 



For most functions yl(W^), this equation does not admit 
solutions expressible in terms of elementary functions. 
For example, consider a function A(VF) which consists of 
more than one term without being linear in W (if it is lin- 
ear then it is integrable in terms of elementary functions 
as already discussed). For some cases of this type it is 
possible to find the solution in terms of elementary func- 
tions by u se of an intrinsic slicing gauge. As discussed 
in section [II B , such a slicing gauge is characterized by 
choosing some simple function of the metric components 
as the independent variable. In the present problem, one 
can reintr oduce the gauge freedom in equations (4.84) 
and ( 4.85| ) by introducing a new independent variable A 

This choice leads to 



such that N = N, 



(taub)a; ^z{Wy 
the decoupled equation 



z{W)-^F{W) , 



(4.86) 



where the dot refers to the new independent variable A. 
Choosing z — F^^^ leads to = A as the independent 
variable (setting the constant of integration to zero), so 
that the slicing gauge is clearly an intrinsic one. Inserting 
W = X into the Hamiltonian constraint and expressing 
this in the new slicing gauge yields 

^ = 2FCxr'[AC\)V + BCX)], (4.87) 

which is easily solved. However, whether or not the so- 
lution can be expressed in terms of elementary functions 
depends on the explicit expressions for A, B and F. 
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To be more exphcit, consider the interior Schwarzschild 
case where the usual Schwarzschild radial coordinate is 
related to the metric scale factors by the intrinsic slicing 
condition r = Ri = R2 = . Using the power variables 

W ^ e"" ^ " 

leads to 



V = §6 3" and the gauge choice x = e^+a" 



H = -IWV + 24[(1 - Kp(o)W^)V + k(p(o) + Pio))W^ , 



st]=F. 



(4.88) 



For the particular case of nonsingular solutions, smooth- 
ness conditions at the center (where W = r ^ 0) (see 
e.g. Q) require that the integration constant occur- 
ring in the expression for the decoupled variable must 
be zero. Reintroducing the gauge freedom and choosing 
z = = [96{W - ^KP(Q)W^)]^^^ leads to the usual 
Schwarzschild gauge r = W = X 

+<P{0) +P{o))r^/[r ~ \>^P(.Q)A ■ (4.89) 

This equation is easily integrated and going back to the 
original metric variables one finds the simple standard ex- 
pression for the interior Schwarzschild solution . The 
general case with a nonzero constant has more compli- 
cated solutions. 



G. Killing tensor symmetries for a subclass of 
2-dimensional models 



By reducing a given problem, either by exploiting sym- 
metries or by specializing to a subcase, one often ends up 
with a reduced system having only a few degrees of free- 
dom. Apart from the trivial case when there is only a 
single degree of freedom left, the simplest reduced sys- 
tems have two degrees of freedom. Many problems can 
be described by a reduced 2-dimensional Hamiltonian of 
the special form 



(4.90) 



where w^v are the standard null variables of Eq. ([1.14| ), 
and 



F{a)., 



(4.91) 



where c is a constant. Although c (if nonzero) can always 
be normalized to unity by a suitable rescaling of a and 
P we choose not to do so here in order to facilitate com- 
parison with the table below. However, the translational 
freedom in (3 will be used to simplify formulas. 

When c ^ the potential form (4.91) corresponds ex- 
actly to the case when the associated Jacobi metric 



admits a homothetic symmetry generated by ^ = d/da 
(timelike homothetic Killing vector (HKV) case) or ^ = 
d/dp (spacelike HKV case) respectively 



£^Jab = '^cJab ■ 



(4.93) 



In the case c — the potential depends only on a sin- 
gle variable and ^ reduces to a Killing vector symmetry. 
It follows from the form of ^ that the above variables 
are adapted to this symmetry. The problem of classi- 
fying t he fun ction F{(3) [or F{a)] for which the Jacobi 
metric (4.92) admits Killing tensor symmetries has been 
analyzed in . As explained in that reference it is suf- 
ficient to consider the timelike HKV case. The spacelike 
HKV case can then easily be obtained by an appropriate 
transformation. 

In this subsection all potentials will be given for 2- 
dimensional models which admit a second rank Killing 
tensor of a given weight under the homothetic symmetry, 
subject to the assumption that the Killing tensor Kab 
is characterized by a homothetic weight 2b through the 
equation £^Kab = 2hcKAB- This includes some cases 
which were not stated explicitly in [^ . 

The classification of potentials admitting such Killing 
tensors depends on two parameters describing properties 
of the Killing tensor. The first parameter is the sign 
of the determinant of the conformal part of the Killing 
tensor, E = sgndet(PAB), where Pab = Kab — \KJab 
and K = K^a- The Killing tensor type is related to E 
according to 



E = 



(null) 

(nonnuU H-J) 
(nonnuU harmonic) 



(4.94) 



correspondin g res pectively to the three cases («), (ii), and 
{Hi) of Eq. (4.10). The second parameter is the homo- 
thetic weight factor b. The cases 6 = 1 or 6 = require 
special treatment compared to 6 7^ 0, 1. With the three 
values of E, this leads to nine different cases altogether. 

We now enumerate the potentials of the form (4.91) 
admitting Killing tensors corresponding to these cases. 
The three cases corresponding to a null Killing tensor 
are collectively given by 



(A) (M1;S = 0): 

U = [two exponential term case I in Table 3.] , 

(B) (6 = 1; E = 0) : [/ = [Cic(w - v) + Caje^'^'" , (4.95) 

where case A also includes the case 6 = = E. When 6 = 
1, the spacelike HKV case is obtained by i nterch anging 
w and V in the corresponding expression in (4.95). 

For nonnuU Killing tensors, the function F(f3) is given 
by one of the following expressions (modulo a translation 
of (3) 



Jab = 2\U\riAB 



(4.92) 



(C) (67^0,1; E = l) : 

Ci cosh" [2c/?/ (s + 2)] -I- C2 sinh"[2c/?/(s -I- 2)] 
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le 



-2c/3(s+2)ls 



[D) (M0,1; S = -l) 

{E) (6= 1; E = 1) : 

[two exponential term case III in Table 3] 

(F) {h=l- S = -l): 

+C2sin[2c(l-A:2)i/2^]}, 

(G) (& = 0;S=1): 
Cilncoth(c/3) +C2, 

{H) (6 = 0;E = -1): 
Ci arctan e^'^'^ + C2 , 



(4.96) 



where s = —26/(6 — 1) and fc (|fc| < 1) are real parame- 
ters while Z? is a complex parameter. In cases (C) and 
{D) the potential can be expressed explicitly as a sum of 
exponential terms if s is an integer. The subcases with 
exac tly tw o exponential terms are given explicitly in sec- 
tion [VH . Case [D] can also be expressed explicitly as 
a real function in the form 

F = cosh'/^[4c/3/(s + 2)]{Ci cos [s arctan e^"'^/ ("+2)] 
-t-C2sin[s arctan e'*"'^/(''+2)]} . (4.97) 

Using multiple angle formulas the trigonometric expres- 
sion inside the curly brackets can be converted to alge- 
braic form provided that s is a rational number m/n. 
However, since this involves solving a polynomial equa- 
tion of degree explicit algebraic expressions can only 
be guaranteed for \n\ < 4. Some of the above Killing ten- 
sor cases admit special Killing vector cases, e.g., setting 
Ci = or C2 = in expression (C) leads to such a case. 



H. Killing tensor cases for Taub potentials with two 
exponential terms 

An important special case of 2-dimensional systems oc- 
curs when the potential is a sum of two exponential terms 



(4.98) 



Extracting all the two exponential term potential cases 
from the various types of the previous section and adding 
the flat case with a null HKV and the nonnuU Killing 
vector case corresponding to c = in Eq. ( 4.9l| ) leads to 
Table 3 for the corresponding parameter values. In the 
case (V) of this table, the type of nonnuU Killing tensor 
depends on the relative sign Z = sgn(CiC2) of the two 
terms in the potential. For Z = 1 one has a nonnuU H-J 
Killing tensor case (S = 1) while for Z = — 1 one has a 
nonnuU harmonic Killing tensor case (S = —1). 

The null (I) and flat (II) cases are easil y treat ed using 
results from the null decoupling section IV E 1 . (Case 



I admits additional nonnuU Killing tensor cases for cer- 
tain parameter values, for such cases one may choose 
nonnuU solution techniques.) Referring to that section, 
the flat case corresponds to Ai — 0, while the null cases 
correspond to Ai ^ 0. In case (III), decoupling can 
be achieved by an appropriate Lorentz transformation 
leading to two Friedmann equations. The case (IV) is 
a nonnuU Killing vector case where decoupling can also 
be accomplished by a Lorentz transformation leading di- 
rectly to a single generalized Friedmann equation. In the 
remaining nonnuU cases (V), one can introduce power 
variables leading to an easily solved problem with a po- 
tential which is a quadratic form in the new variables. 
An example of such a case has been dealt with in p3| . 



V. INVARIANT SUBMANIFOLDS AND HOW TO 
OBTAIN THEM 

Solving the Einstein field equations in general seems to 
be impossible, particularly in view of recent results that 
the only generalized local symmetries of these equations 
arc due to scale invariance and the diffeomorphism group 
[Q, and these symmetries are insufficient to lead to a 
general solution. To find special solutions one imposes 
space-time symmetries and/or other restrictions on the 
dependent variables so that one obtains a more tractable 
consistent subsystem of differential equations. In other 
words one tries to find "invariant submanifolds" of the 
original system of field equations. Even imposing enough 
space-time symmetries to reduce the field equations to or- 
dinary differential equations as one does to obtain the HH 
models still does not lead to such tractable subsystems 
in general. One must impose further conditions to be 
able to actually find exact solutions. There is no general 
systematic method of discovering invariant submanifolds. 
It is here that creativity and imagination and even plain 
luck play a role in rooting out these hidden structures. 
There are many particular ways in which invariant sub- 
manifolds have been found, but few of these successes 
involve a systematic method. Many methods require an 
arbitrary function like an unspecified equation of state 
or an unspecified scalar field potential to produce solu- 
tions. In this brief section one systematic method will 
be presented which does not rely on the existence of ar- 
bitrary functions and is relevant to many though not all 
of the known invariant submanifolds. In particular for 
Hamiltonian problems this method also yields the class 
of exact power law (EPL) solutions. EPL solutions have 
been studied in p5| , p6| . 

Hamiltonians which are reducible to the following form 
play a crucial role in the discussion of HH models 



H = 



(5.1) 



where x and U are analogous to the previous slicing gauge 
function x and the Taub potential. Choosing x = 1 leads 
to the equations 
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ij'' = -r]'"'dU/dy'' . (5.2) 

If dU/dy^ = holds for some value y'^^^ of a par- 
ticular coordinate y^ independent of the values of the 
remaining coordinates, and if this condition is compati- 
ble with the Hamiltonian constraint, then y'^ = de- 
scribes an invariant submanifold. The equations for the 
remaining variables on this submanifold are given by the 
above Hamiltonian after having inserted the conditions 
— 2/fo) ^^^"^ ^ (^'-"' ^^^^ coordinate alone). Invari- 
ant submanifolds within invariant submanifolds are also 
possible. Reduction down to one dimension automati- 
cally leads to a solution since the Hamiltonian constraint 
only involves a single variable, thus leading to a quadra- 
ture. 

Lorentz transformations of the variables are sym- 
metry transformations of the Minkowski metric appear- 
ing in the expression for the kinetic energy function. 
They play a crucial role in finding many invariant sub- 
manifolds. All HH models of the previous section have 
Hamiltonians or reduced Hamiltonians with Taub po- 
tentials which can be put into the following form by a 
Lorentz transformation from the variables to new ones 

i 

If dFi I dp^ = holds for all values of i for some particu- 
lar value /3^') of a particular coordinate independent 
of the values of the remaining coordinates, and if this 
condition is compatible with the Hamiltonian constraint, 
then = P^Q.^ describes an invariant submanifold. In 
many HH cases the index value and some definite value 
P can be interchanged in this discussion, but the 
Hamiltonian constraint seems to prevent the existence of 
invariant submanifolds of this type. 1-dimensional invari- 
ant submanifolds with one exponential term lead directly 
to EPL solutions. 



(Killing vectors and/or homothetic Killing vectors) or 
additional continuous intrinsic symmetries (Killing vec- 
tors). The latter are symmetries of the intrinsic geom- 
etry of the individual homogeneous hypersurfaces which 
are not necessarily spacetime symmetries. As in section 
2, the diagonal and nondiagonal models are treated sepa- 
rately, but the diagonal models are collected according to 
the dimension of the intrinsic symmetry group. Unless 
otherwise stated, the only perfect fluid solutions being 
considered here are those for which p — — l)p. 

A. Diagonal models 

The possible dimensions of the intrinsic symmetry 
group of the geometry of the HH hypersurfaces are 6, 4, 
and 3. Beginning with dimension 6, models are consid- 
ered with only /3'^-dependent sources and possible scalar 
fields. This class of models includes the Bianchi type 
I and V models and the SH constant spatial curvature 
type IX models (there are no static models of this latter 
type). The SH models belonging to this class are intrin- 
sically isotropic. 

Next diagonal models with a 4-dimcnsional intrinsic 
symmetry group and with only /J*', /^"'"-dependent sources 
are treated. These models are all intrinsically LRS and 
include the Bianchi type I, II, and V models, the LRS 
Bianchi type III, VIII, and IX models, the SH Kantowski- 
Sachs models, and the static spherically symmetric mod- 
els. Note that apart from the SH Bianchi type IX FRW 
perfect fluid solutions and the SH LRS Bianchi type VIII 
and IX stiff perfect fluid solutions, there arc no other 
known exact perfect fluid solutions for these two Bianchi 
types. 

Finally the diagonal SH Bianchi type VI vacuum and 
perfect fluid models are considered. 

1. Models with a 6-dimensional intrinsic symmetry group 



VI. PROBLEMS LEADING TO EXACT 
SOLUTIONS AND HOW TO SOLVE THEM 

This section will survey the cases which lead to exact 
solutions. The method of solution which works in each 
case will be specified by referring to previous sections, 
without going through the mechanical details of obtain- 
ing and presenting the solution explicitly. In fact, as has 
been shown, there are often several ways one can solve a 
given problem and hence more than one representation 
of the solution exists. Specific examples of how to use the 
methods of this article to produce the actual spacetime 
metrics which correspond to these solutions are given by 
Uggla and Uggla and Rosquist |37|] . 

Except for a few special class B cases and Bianchi type 
VIo, all exact solutions arise from spacetimes which ad- 
mit either additional continuous spacetime symmetries 



a. Sources not including a scalar field. The Hamilto- 
nian can be written as ||38|| 

H=\x{~p''^+P+^+p-^) 

+x-^[~72ke^P'' + U(,,^^,){P°)], (6.1) 

where C/(taubs) is the source potential expressed in the 
Taub slicing gauge. The variables are equal to zero 
for the type IX models while P'^ is equal to zero for the 
type V models. One can choose 1-fornis so that the pa- 
rameter k has the values for Bianchi type I, 1 for type 
IX (correponding to the choice n^^^ — rS^'^ — n^^^ — 2), 
and —1 for type V (corresponding to the choice a — 1). 

Letting x only depend on results in the equations 
of motion {xP^)' — 0, which lead to constant values of 
the momenta p± = xP^ . Thus are determined by the 
equations 
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/3± - X- V± • (6.2) 
One may solve for using the Hamiltonian constraint 



(/3°)2 = .T-2[I]2-72fce4'3" + [/(^ 



aubs) 



(6.3) 



where Y? = p^"^ + pJ^ . This equation immediately gives 
a quadrature for The most interesting case is when 
U{taxLh&) is a sum of exponential terms and this problem 
reduces to the the generalized Friedmann problem. Ref- 
erences to some of the literature on the most notable SH 
solutions with a 6-dimensional intrinsic symmetry group 
are given in Table 4. The SH vacuum type I solution 
is usually associated with Kasner who found the corre- 
sponding static solution |^ . The isotropic vacuum type 
V solution is just the Milne form of Minkowski spacetime. 
Useful references for FRW and FRW-A models are Har- 
rison go), Vajk [|l[, Anderson Q and Misnei% Thorne 
and Wheeler The book by Kramer et al [Q is also 
useful in this context as well as for further references on 
models in Table 4 with symmetry groups of dimension 3 
and 4. 

h. Sources including a scalar field. The Hamiltonian 



-2Ax~^[ike^l^° + eKeS^V(3c)(/3^)] ■ (6.4) 



Solvable cases. For most scalar potentials this is 
not a solvable problem. However, if V(sc) is the sum of 
exponential terms, then some solutions do exist. An in- 
teresting example is the case of scalar field models with 



-2C/3T 



For 



a single exponential potential V^g 
such models with = 0, which includes the isotropic 
models, the Taub potential is given by 

f/(taub) = -24[3fce2""+2^' + eKe(3-'=)™+(3+=)^] , (6.5) 

where w = (3^ + [3^ and v = [3° - 13^. If fc = this is 
a simply solvable one-exponential-term problem. When 
fc 7^ there are two solvable cases. The first case c = 
1 correspo nds to the Ai — Q flat null decoupling case 
of section IV E 1 . The second case c — 2 corresponds 
to the nonflat null Killing tensor case The Jacobi 

metric of the case fc = and (3^ — with an arbitrary 
scalar potential, V(sc~i , admits a timclike HKV. Therefore 



the Killing tensor cases {A) through [H) of section [V G 
apply and lead to exact solutions. The particular case 
(C) with s = 2 leads to the solutions found in [^ , p7| . 
However, one can easily produce many other solutions of 
comparable physical interest. 

Invariant sub manifolds. If V(^sc) ha^s relative ex- 
trema, then one has an invariant submanifold corre- 
sponding to the corresponding flxed value of the scalar 



field. The resulting problem yields a generalized Fried- 
mann equation where the scalar potential reduces to an 
effective cosmological constant. 

There are other more interesting invariant subman- 

j58|j5|@|ll. 



ifolds obtained by a different method 
These correspond to exact solutions describing inflation- 
ary models in cosmology as well as static domain walls 
in an astrophysical context. 



2. Models with a ^-dimensional intrinsic symmetry group 

For the family of intrinsically LRS class A models, 
which can be chosen to satisfy n^^^ = rS^^ , it is convenient 
to introduce the notation a = n'^^-'n^^-'. One must set 
n^-^) = to obtain the remaining LRS models, for which 
the curvature parameter a continues to have its previous 
meaning. The sources considered in this subsection may 
include a cosmological constant, electro- magnetic fields 
and perfect fluids. 

The Hamiltonian for this family of spacetimes is given 

by 



[A 5 = 0,+,-], (6.6) 

where /3~ must vanish except for the Bianchi types I and 
II where it is a cyclic variable, provided that x is assumed 
to be independent of this variable. The system associated 
with this Hamiltonian provides several interesting exam- 
ples of null and/or nonnuU Killing tensor cases which in 
turn give rise to many exact solutions. 

a. Vacuum, A, EM field, SH stiff perfect fluid. 

The case A = 0. The Hamiltonian for this case can 
be nicely expressed in terms of a new pair of conformally 
inertial coordinates obtained by the Lorentz transforma- 
tion § 

(^°,^+) = 3-1/2(2/3" -/?+, -/30 + 2/?+) , 

(/3°,/3+) = 3-1/2(2^°+^+,-^° -2^+), (6.7) 

in terms of which the Hamiltonian takes the following 
form 



H ^ ix(-^" ' +/?+' + /?"') + 2Ax-^[\n^^^^ 2^-473^+ 

(6.8) 



It is easy to see that this Hamiltonian is a nonnuU H- 
J Killing tensor case. Furthermore in the Taub slicing 
gauge X = I, one has a completely decoupled Hamilto- 
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(taub) 



- + Hj, 



13" 



H0- 



/3+ 



(stiff) I 



l„(3)2g-4V3?+ 



(stiff) 



24kp(o) -B(stiff) (e -1) 



2 - 

ee e 



/3+ 



(6.9) 



The full problem therefore consists of three generalized 
Friedmann problems interpretable as 1-dimensional prob- 
lems with exponential potentials and constant energies, 
restricted only by the constraint 



Ero + Et,+ 



Ep- 



E, 



(stiff) 



0. 



(6.10) 



Each of the 1-dimensional motion problems is governed 
by a generalized Friedmann equation, with a potential 
for the variables (when nonzero), l3^, and hav- 
ing in general one, two and three exponential terms re- 
spectively, the latter occurring as an "equally spaced" 
exponential coefficient case (see section |l|), all of which 
are equivalent to generalized Friedmann problems with a 
quadratic potential. 

Apart from the trivial cyclic variable (3~ which is 
present for Bianchi types I and II, the natural variables 
which lead to quadratic potentials for the other two de- 
grees of freedom are 



U- 



/3+ 



g2v^+ ^ ^(3) ^ 
3/3+ „(3) ^ . 



(6.11) 



The case /?- = = JJ, 



(fluid)- 



These models may be 



re-examined as an example of a 2-dim ensional null Killing 
tensor case treated in section IV E 1 , Letting 



w^l3" + l3+ , v^p^ -13+ , 
the Taub time gauge potential takes the form 



U, 



(taub) 



_gg2g-.-+3. 



g^g3(^.+«)] 



(6.12) 



(6.13) 



For the vacuu m ca se this expression can be identified 
with equation ( |4.47| ), with Co = 3 and Aq = 2{n'-^yf/3, 
allowing w to decouple. This permits a nonzero cosmo- 
logical constant term since only a single null variable is 
required to decouple, in contrast with the previous non- 
null discussion where the cosmological constant had to 
be zero. 

The natural power variables and slicing gauge function 
are 



V = e' 
leading to 



3u 



W 



X = imVW , 



(6.14) 



H 



8.rl„(3)2 



eAI^m-i] 



(6.15) 



Since Oq = — 2 7^ 0, the three choices m = —2, ±1 lead 
to to elementary function solutions for W . The choice 
TO = 1, which makes the e term a constant, was first 
introduced in the SH context by Misner and Taub 0. 
The choice m = — 1, which makes the term a constant, 
was first introduced by Brill in the SH case. A third 
new slicing gauge arises for the choice to = — 2 which 
makes the n'^^^ term proportional to V . Note that — 
for the spherically symmetric models. Thus the choice 
TO = 1 leads directly to the standard expression for the 
Reissner-Nordstrom solution with cosmological constant. 
If one also sets A = then one has a nonnuU-nuU case 



which may be solved as in section [VE3 . 

If there are several terms equal to zero there are even 
more slicing gauges and dependent variables one can 
choose to solve the problem. As an example, consider 
the LRS Bianchi type II, III and the KS vacuum mod- 
els which correspond to a 2-dimensional problem with a 
Taub potential consisting of a single exponential term. 
This term is a nonnuU exponential, and the correspond- 
ing pro bl em is ea sily solved using the methods of sections 



IVEl or IVE2 



References to some of the literature on the more promi- 
nent solutions are given in Table 5. Apart from the so- 
lutions indicated in this table, it's worth noting that the 
general LRS solution with an electromagnetic field and a 
cosmological constant have been given by Cahen and De- 
frise [Q. A useful reference and guide to the literature 
on solutions with electromagnetic fields is the work by 
MacCallum |Q. MacCallum, together with Siklos, has 
also made a thorough investigation of HH vacuum mod- 
els with a cosmological constant . For a discussion on 
LRS models see |Q. 

h. Static perfect fluids. The most interesting static 
models are the astrophysically relevant spherically sym- 
metric ones. The Bianchi type I models are also of some 
interest as cylindrically or plane symmetric (/3~ — 0) 
models. Other static Bianchi models do not seem to be 
particularly interesting physically and will not be consid- 
ered here. 

Spherically symmetric models. For the astro- 
physical spherically symmetric models, various equations 
of state have been considered. 

The case p = (7 — l)p : For the usual equation of state 
with 1 < 7 < 2 one has the Taub potential 

C/(taub) = 24[e4^"-2^" + «p(o)e(«-'')^°+2"^"] , (6.16) 

where 77 — 7/(7 — 1). 
Making the boost 

/30=r(^" + z;^+), 

/3+ =r(w;5°+^+), where r = (1 - t;2)-i/2 , (6.17) 
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with the value u = ^{t] - 2)/{r] + 1) = i(2 - 7)7(27 - 1) 
of the boost parameter, leads to 



A 



(taub) 



3r(>7+2) 



B 



er 



C 



2(>)+l) 



(6.18) 



This potential is of the same form as Eq. ( |5.3| ) with a 
single term and has a nonzero minimum value. A 1- 
dimensional invariant submanifold corresponds to this 
minimum value, leading to a generalized Friedmann 
equation with one potential term easily solved using the 
methods of section ID. This solution is a special case 
of solutions found by Tolman The solution with 

7 = 4/3 has also been found by Klein [ p3| . 

Th e case p — /9(o) : This case has been treated in section 
[VF. The nonsingular interior solution was first found 
Solutions which have a singularity 



by Schwarzschild [[ 
are also null cases and have been investigated by Volkoff 
Q and Wyman [|5|. 

The case of an unspecified equation of state: The Taub 
potential is 



(taub) 



24[g4/3«-2/3+ ^ ^e^0"p(^po _ 2^+)] _ (6.19) 



Making the boost ( |6.7D leads to 



24[e2v^" + ^e4v^V(^+)] , (6.20) 



(taub) 



where / = e^^'^" p(l3~^). If / has a minimum for some 
value of /3+ then one obtains an invariant submanifold as 
discussed in section ^ Unfortunately such a minimum 
leads to an unphysical equation of state p = — p/3. 

However, this problem has the same form as the scalar 
field problem with an unspecified scalar potential men- 
tioned above and dealt with in [Q. Thus the same 
method can be applied to the present problem and will 
produce invariant submanifolds and corresponding ex- 
act solutions. Alternatively one can specify f{f3~^) or 
p(/3+) to be some function so that one obtains a problem 
for which one might find an invariant submanifold or a 
Killing tensor and thus exact solutions. Once a solution 
is found the equation of state can be derived. Unfortu- 
nately, the general solution to the Killing tensor problem 
is not available at present. However, all of the solvable 
cases found in the literature can be recovered by a cer- 
tain ansatz for a conformal transformation relating the 
standard null variables {w, v) to a set of null variables 
{W, V) which are adapted to the symmetry Q . 

The starting point is to write down the Jacobi metric 
in standard null variables w — f3^ + /3+, v — (f' ~ I3~ 
leading to (modulo a constant factor) 

ds/ = -2 {e"'+-''^' + Ke3(-+'')p(^i£±32i)| dzwdw . (6.21) 

The ansatz we use for the conformal transformation is 



where r and s are constants to be determi ned. Applying 
this transformation to the Jacobi metric (3.21) yields 



-2G{W, V)dWdV , 
G = rs [W-^V^"-^ + KW'^''-^V^''-^h{Y)\ 



(6.23) 



where h{Y) = p{\ogY) and y = e'^' = VF^'^V^^^/^ {j^c&W 
that (3^ ~ (3^ — 2/3+). We next look for conditions on r, 
s and h{Y) which make the variables W and V symme- 
try adapted with respect to a Killing vector or a Killing 
t enso r. We do this by in sertin g the expression for G in 
( |6.23 ) into the equations (4.10). Analysis of the resulting 
set of equations leads to the solvable cases given in Table 
6 (for how one exphcitly solves these cases see Q). 

The cases for which the equation of state is of physical 
interest are the null case with s = 2/3 (Schwarzschild's 
interior solution), the Hamilton- Jacobi case with r = 1, 
s = 1/3, (the Killing vector case, where a — b, corre- 
sponds to Buchdahl's generalized polytropic solution of 
index five [^) and finally the Hamilton- Jacobi case with 
r = 2, s = 2/3 (setting c+ = gives Buchdahl's gener- 
alized polytrope of index one [^,^, while the general 
c_ 7^ case was recently given by Simon [pO|). 

Bianchi type I models. Again several equations of 
state are of interest. 

The case p = (7 — l)p: These models correspond to a 
problem with a single exponential potential term 



C^(taub) = 24kP(o) 



(6.24) 



which is nonnuU for 1 < 7 < 2 and null for 7 = 2 since 
7] — 7/(7 ~ 1). These two types of problems were dealt 



with in section [V E 1 and IV E 2 



The case p — p(o) + (^7 ^ ^)P- The Taub potential is 

f/(taub) = 24.[(p(o) - P(o)e«^"] . 

(6.25) 

For rj — 6 (or 7 = 6/5) it is easily seen that one has a 
nonnuU Killing tensor case, thus easily solved with the 
methods of section IVE 2| . The solution was first found 
by Evans 

c. Spatially homogeneous Bianchi type II nonstiff per- 
fect fluid models. These models have the Hamiltonian 

+2Ax-\\e-^^^^^-P''^ + «;p(o)e3(2-7)/3''] , (6.26) 



where [A, B = 0, +, -] . Making the boost ( |6.17D with 

« = i(37-2) 
leads to the Taub potential 



(6.27) 



(6.22) 
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C/(taub)=6e3r(2-7)/5''[e-3r(6-.)/3+/2 



-4KP(o)e 



er(2-7)(37-2)/3+/8l 



(6.28) 



The /^"''-dependent factor in this potential has a min- 
imum. As discussed in section |V this corresponds to 
a 2-dimensional invariant submanifold, leading to a 2- 
dimensional problem with on e nonn uU ex ponent ial term, 
easily solved using sections IV E S and [V E 1 . When 
^ this yields the Collins solution |92[, while for 
= it gives an LRS EPL solution [||,g3|]94). 

d. Spatially homogeneous KS and Bianchi type III non- 
stiff p erfect fluid models. Referring to sections 11 C 3 and 
11 C % the orthogonal perfect fluid models have a Hamil- 
tonian of the form 



(6.29) 



Transforming to the null variables w — (3^ + /?+ and 
= Jf' — allows the potential to be identified with 



(4.98) of the two term case of section [VH, with the 
following correspondence between the parameters 

Ci = 24a, C2 = 24«p(o), 

Pi = l,gi = 3, P2 = 92 = 3(2 - 7)/2 . (6.30) 



Solvable cases occur for the following parameter values: 

The fiat null = case. Condition (//) in Table 
3 corresponds to the radiation value 7 — 4/3. 

The null Ai 7^ case. Condition (/) yields physical 
solutions for the dust value 7 = 1 and the value 7 = |. 



The above solutions 
|l|,||,|l9i9||3. 



can be found 



3. Spatially homogeneous Bianchi type VI models 

a. Vacuum models. All known type VI vacuum mod- 
els are Taub symmetric {f3~ — 0) and correspond to a 2- 
dimensional problem with a Taub potential which is the 
single exponential term in equation (2.2C), which is non- 
null except for Bianchi type VIq where it is null. Th ese 
are eas ily solved using the methods of sections IV E 1 and 



IVE2| . These solutions were first found by [||j9l]j9|,0 



b. S olv able pe rfect fluid models. Referring to sections 
II C 3| and II C 2 , the orthogonal perfect fiuid models have 
a Hamiltonian of the form 

i7=ia;(-/30 2+/3x2) 

+2Ax-'[c-~^e^^^''--^^^^ + «p(o)e3(2-7)/5°] , (6.31) 
where c"^ = + Sa^. 



Transforming to the null variables w = (3^ + f3^ and 
V ^_Jf' — allows the potential to be iden tified with 



( 4.98 ) of the two term case of section IV H , with the 
following correspondence between the parameters 

Ci = 24c-2 , C2 = 24kp(o) , pi 2(1 - eg) , 
qi = 2(1 + cq), P2 =92 -3(2-7)/2. (6.32) 

The physical cases correspond to 

0<cg<l, 1<7<2. (6.33) 

Solvable cases occur for the following parameter values: 



The flat null Ai = case. Condition (/) in Table 
3 yields cq = |(37 — 2). The resulting solution was first 
found by Collins and has been presented by Wain- 
wrigh t [ p5[ corresponding to the form given in equation 
(4.55) with m = §. 



The null A., ^ case. Conditions (/) of Table 3 
with ao 7^ yield the conditions cq = ±^(4 — 37) or 
eg ~ i(37-|-2), respectively, and solutions found by Uggla 
p3| and Uggla and Rosquist |]37|| . 

The nonnuU case. There also exist two examples 
of the nonstiff perfect fluid solutions corresponding to 
nonnuU cases in which the potential may be reduced to a 
quadratic expression in the two natural power variables 
(case (V) in Table 3. These correspond to the values 
eg = I and 7 = | , leading to the solution given by Uggla 
p3[ . These exact solutions together with the above null 
type VI^ ones were found using the methods developed 
in the present article and are the first new orthogonal SH 
non-EPL perfect fluid solutions found in several decades. 
The type Vlh {h ^ 0) stiff perfect fluid solution can be 
obtained by making a boost that leads to a generalized 
Friedmann problem (the type Vlh=o stiff perfect fluid 
case is contained in the flat Ai = case discussed above) . 

Invariant sub man ifold perfect fluid models. 

Making the boost ( |6.17| ) with value v = (37- 2)/(4cg) < 
1 leads to the Taub potential 



C/(taub) = 24e3r(2-7)/3"[c-2e[37-2-4(c,)^]/3X/(c9) 



-HkP(o) 



,3r(2-7)(37-2)/3V(4cg)i 



(6.34) 



For the type Vlh models, the -dependent factor of 
this potential has a mininmm provided that (37 — 2)/4 < 
{cq)^, yielding an EPL solution. This solution was dis- 
cussed by Hsu and Wainwright ||5l[| . 

Bianchi type VIo mod els. For nonstiff perfect 
fluid models the boost ( |6.17 ) in the /?+-direction with 
the value v = — 2) leads to the Taub potential 
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C/(taub) = 6e'^^'-''^^°[e^^^h^'+Anp^,^e^~^' 



A = 3r(2-7), B 



[r(2-7)(37- 2). (6.35) 



The /S"*", /3~ -dependent factor has a minimum, yielding a 
1-dimensional invariant submanifold which corresponds 
to the same type VIq EPL solution just discussed. How- 
ever, these new dependent variables are useful for a qual- 
itative discussion of the dynamics for this class of models 



B. Nondiagonal models 

1. Stationary cylindrically symmetric models 

The stationary cylindrically symmetric vacuum models 
have the reduced Hamiltonian ( ^.50| ) which corresponds 
to the single nonnuU exponential potential te rm cas e and 
is therefore solvable as discussed in section IVE 2 . For 
an explicit representation of this solution sec |4|. 



2. Spatially homogeneous Bianchi type Vl-i/g models 



a. Vacuum models. The boost (S.17) with the value 
V = 2/(3V3) transforms the Taub potential from ( 2.56 ) 
to 



„8r/3V3r3„ 2 4\/3r/3+ 



32e 



-ioy3r/3+/9i 



(6.36) 



The /^"'"-dependent factor has a minimum, which corre- 
sponds to a 1-dimensional invariant submanifold, leading 
to an EPL solution [0195 



b. Perfect fluid models. The same boost as in the 
previous vacuu m cas e for j — 10/9 transforms the Taub 
potential from (2.5(:) to 



U, 



taub) 



^24^io/9e8r/3+/3i 



(6.37) 



The /J^-dependent factor has a minimum, which corre- 
sponds to a 1-dimensional invariant submanifold, leading 
to an EPL solution . No known exact solutions exist 
for other values of 7. 



3. Spatially homogeneous class A models belonging to the 
symmetric case 

a. Bianchi type II perfect fluid models. Here we are 
going to show how one can obtain tilted type H EPL solu- 
tions by finding invariant submanifolds without explicitly 
knowing the function Y occuring in the fluid potential. 
Tilted models are quite complicated and therefore the 



manipulations become rather cumbersome. However, by 
using comput or alg e bra th ey ca n be done. 

Equations ( ^.16 ), ( 2.67 ) and ( 2.73 ) yield the total po- 
tential 



C^(taub) — C^(c) 



C^(G) + fluid) 



(6.38) 



+ Se4(/3»+/3+ + v^/3-) ^ Cg3(2-7)/3V-^/2(^r " 7 + 1)] , 

where B and C are constants defined by 



C = K-i(n(i))2r(2-^)(z;3)- 



(6.39) 



It can be simplified by first performing a boost with ve- 
locity V — (4— 37)/2 (excluding the stiff fluid case 7 = 2) 
in the (3^^ direction 

/3" = r,(^° + t-/3+) , /3+ = r„(t;/3° + P+) , p- =p-, 

(6.40) 

where 

r. = (1 - i;2)-i/2 ^ 2[3(37 - 2)(2 - j)]-'/^ . (6.41) 



Then the relation (2.6£) simplifies to 
where ^ is a constant given by 



(6.42) 



(6.43) 



thus leading to y = Y{/3+) since FV-^ - y + 1 = 0. 
Expressed in these variables the potential takes the form 

C/(taub) = 24.^2^2(^^)2(„(l))-2[g-4V^r (6.44) 
_^(^g(3/2)(2-7)r„[2/3«-(37-4),3+]y-7/2(^y - 7 + 1)] . 

A further boost with velocity w = — (\/3/4)r^(2 — 7) = 
-V^~^/[2V37- 2J in the direction 



/3° = r„(/30+u;/3'), /3+=/3+ 



where 



r^ = n-,„2)-i/2= ^ ^^^ 
^ ' V137 - 10 



(6.45) 



(6.46) 



leads to 

[/(taub) = 24«2^2(.3)2(n(i))-2e3(2--)r"'^"$(r,r ) 



(6.47) 
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where 

+ (7e(3/2)(2-7)(4-37)^+-4\/3r„(l+a)/5-y-7/2('^y - 7 + 1) . 

and a — — (r^,)^^. This form of the potential displays ex- 
plicitly its homothetic character in the sense that d/dfr 
is a homothetic vector of the associated Jacobi geometry. 

To locate possible extremal values of the function $ it 
is convenient to express it as 

+CA-^Z^+'^Y''{Y - lf{-iY - 7 + 1) , (6.49) 

where we have introduced the variable Z = e^^^/'"™)'^ 
and the constants 

57-4 , 4-87 

0. = —-ri-: , = 



k ^ a 



2(37-2)' 2(87-2)' 

4(37-2) 



-(r„ 



137 - 10 



(6.50) 



The potential can be further simplified by making a suit- 
able translation in P~ or equivalently rescaling Z by 
Z = SZ where S = jsh+\l2Q-\i2 ^ rpj^ig gi^^g s^^^^^^ 

$ = J'^'*!) where 



$ = Z*^' + D{Y)Z^^'^ + E{Y)Z 



7k+l 



where 



D{Y)^Y'^''+\Y -lf^+^ , 

E{Y) ^ 2r"(r - i)\Y - 1 + 7-I) , 



(6.51) 



(6.52) 



The conditions Z > and Y > I can only be fulfilled if 
7 > 10/7. 

We conclude that there exists an exact power law so- 
^gl|i^i|On with values of corresponding to the minimum 
values of Y and Z through the relations 

g(4/r„)/3+ ^ ^-ip ^ ^-iyi-7(y _ 1) ^ 

g(4V3/r„)r ^ ^6+1/2^-1/2^ (g 

This is exactly Hewitt's EPL solution for the equation of 
state parameter in the range 10/7 < 7 < 2 p9|l. (Note 
that the solution with 7 = 10/7 occurring in P9|| is not 
a tilted one.) Since the Jacobi metric admits a timelike 
HKV, one can apply the methods developed in ||9^ to 
obtain a complete picture of the behaviour at late times 
of these models forlO/7<7<2. 

b. Bianchi type VIq stijf perfect fluid models. For stiff 
perfect fiuids and the usual values of the structure con- 



To find possible extremal points of $ we calculate the 
derivatives 

$ y = Z'^+i [D'{Y)Z + E'{Y)] , 

^ ^ = Z'^-^ [k + D{Y)Z^ + E{Y)Z] , (6.53) 

where 

D'{Y) = Y'^''{Y - lf\2{a + 6 + l)r - 2a - 1] , (6.54) 
E'{Y) = 2r"-iy''-^{(a + b + l)Y^ 

+ [-2a + 5 - 1 + (a - b)^-^]Y + a(l - 7-^)} . 

Equating the first of these expressions to zero while not- 
ing that D'{Y) > gives Z = - E' {Y) / D' (Y) . Inserting 
this result in the second equation yields an equation for 
Y having a single solution given by 



Y -I 



(37-4)(77- 10) 



leading to 

Z ^Y-''{Y -\)-^ 



27(177-18) ' 
47(177-18) 



(II7- 10)(77- 10) 



(6.55) 



(6.56) 



stants ri 
form 



(1) - 



i'-^-' = 1, the Taub potential takes the 



f/(taub)=6(4 + «:z;3')e^(''°+''")+24«:£2 



(6.58) 



Transforming to null variables this is seen to corr espond 



to the easily solved flat null case co = of section IV E 1 



This solution was first found by Wainwright et al llOi 



C. Summary of the spatially homogeneous perfect 
fluid models 

All of the exact SH perfect fluid solutions arising from 
Hamiltonian models are now collected in Table 7. For 
the sake of completeness the remaining known "non- 
Hamiltonian" exact tilted perfect fluid solutions are listed 
in Table 8. Whether or not some of the models cor- 
responding to these cases admit a Hamiltonian of the 
"standard" form H ~ T + U is not clear. We have not 
been able to find such a formulation for any of the mod- 
els using the present approach. Perhaps the "comoving" 
approach of MacCallum and Taub |101| could be used to 
produce Hamiltonians in at least the tilted type V and 
type VIo cases. 



VII. BEYOND 4-DIMENSIONAL GENERAL 
RELATIVITY 

The techniques presented in this article are valuable 
even outside the realm of 4-dimensional general relativity. 
The only requirement needed to apply them in this larger 
context is a Hamiltonian theory for which the Hamilto- 
nian has a quadratic form kinetic energy function. 
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A. Higher dimensionsional theories 

With the emergence of unified field theories con- 
structed using higher-dimensional spacetimes, the door 
opened to the search for higher-dimensional analogs of 
many of the 4-dimensional HH solutions of Einstein's 
equations or related field equations. These occur as solu- 
tions of Einstein's equations in higher dimensions and its 
numerous generalizations — Kaluza Klein and supersym- 
metric variations of Einstein's equations, with coupling 
to various matter sources — both in the context of static 
HH spacetimes and SH cosmological spacetimes. 

Although details change, the same general picture ap- 
plies and again one sees in the literature on this topic the 
same kinds of similarities that characterize those HH so- 
lutions listed by Kramer et al in the 4-dimensional case. 
The simplest models for which solutions can be found are 
again diagonal, most characterized by the natural gener- 
alizations of the 4-dimensional Hamiltonians, with the 
same Lorentz structure of the kinetic part playing a key 
role in the properties of the dynamics. One may take 
the mathematical discussion of the present article and 
apply it with slight modifications directly to the higher- 
dimensional case [110|. 



B. Nonminimally coupled scalar fields 

Nonminimally coupled scalar field models are de- 
scribed by an action of the form 



(7.1) 



where ^^^g = dei{gap)- If one is so inclined, one may 
also add other matter fields to this action, and below 
we will add a perfect fluid contribution. The function 
B{(j)) can be set equal to 1 without loss of generality 
by using the freedom to redefine the scalar field by a 
transformation (j) — f{(f)) together with an accompanying 
redefinition of A. However, it is convenient to retain a 
general expression B{(f)) for easier comparison with the 
literature. The coupling between gravity and matter is 
said to be minimal if A{(j)) is constant. One can then set 
A = 1, B = 1. All other couplings are referred to as 
nonminimal. Often used nonminimal couplings include 
conformal coupling characterized by A{(f>) = 1 — ^0^, 
i3((/)) = 1, and Brans-Dicke couplingfor which A(0) = 
K(j), B{(j>) = (^/(p, and V{cj)) = 0. Some typical scalar 
potentials are 



y(0) = 



A/k 

2 —2 



cosmological constant, 
mass term, 
Ac/)** theory, 
exponential potential. 



(7.2) 



where m, A and are the physical parameters character- 
izing the various cases. 



Varying the action (7J) with respect to the metric and 
the scalar field yields 

5S 



Sgap - |v/^(«"'^(0)Ga/3 + AJ-i(nA(0)).9„^ 
-K~^A{<^),^p - B(0)0,„<^,^ -I- \\B{^){y<i>f -f VmSo^p) ■ 



^ - v/^[B('^) □ + (2k)-M'(0) Wi? 
-V\^)-\B'm^^f\, (7.3) 
where (V/)^ = .g"^'/,^/,^ and □ / - 5"''/;a/3. 

1. Diagonal spatially homogeneous scalar field models 



We will assume that the scalar field is SH and therefore 
a function of t only. Adding an orthogonal perfect fluid 
with an equation of state p = (7 — l)p to the action (7.1) 
for the diagonal SH models, leads to the Hamiltonian 



SS 



Sg 



2KNg^/^n^nl'Tf^^i)^p (7.4) 



+N 



13/3" 



A{cp)r,ABp^P'' - ^'(0)0/3° + ^B 



V*A((I)) + 2^6^^ V{(j)) -f 2k/9(o)C 



6 

3(7-l)/3° 



where T(pf )ct^ is the energy-mo ment um tensor of the per- 
fect fluid. Here V* is given by ( l.lt ) for the class A mod- 
els for which A,B = 0,-t-,-, by V* = 2c-^e^°-^'"''^'' , 
A, B = 0, X for the type V and VI/i models, by V* — 
-2CTe'^''-2/3+^ ^ 0,-1- for the MT models, and by 

V* = —6k for the FRW models, which are described by 
the line element 



1 — kr-^ 



r\de^ 



(7.5) 



with k taking the standard values 1, 0, —1 for the closed, 
flat, and open models respectively and for which Ini? = 
/3° , = 0. (The isotropic case can of course be ob- 
tained from the other models. However, since this case 
is by far the most discussed when it comes to nonmin- 
imally coupled models, it is given explicitly here.) The 
kinetic energy of a minimally coupled model can be trans- 
formed to a manifestly conformally flat form by the choice 
B{(j)) = constant corresponding to a redefinition of </>. 
For nonminima lly c oupled models we make the field re- 
definitions (c/. pill) 



/30 = /?" + lnAi/2, /3^ = /3^, (7.6) 
- 12-1/2 j ^Z[A'{4>)f + 2kA{4>)B{4>) , 



30 



To appear in Phys. Rev. D 



Preprint, March 1995, 



gr-qc/9503061 



where P takes appropiate values, e.g., ± for the class A 
models. This transformation leads to a manifestly con- 
formally flat kinetic energy 

T = GN-'e'P" [i(/3t)]-i/2 (r^^J^pB + 'pi 2) 

= ^X(„™) (77^5/3^/3'' +/3^'), (7.7) 

where i(/3t) = A{(j>) and a;(„„) = UN-^e^^" [A{p'')]-'^/^ 
— A{(f>)x. Note also that the redefined scalar field is 
consistent with the definition made in section |l| so that 
= in the minimal coupling limit A ^ 1, _B ^ 1. 
With the above variable choice the potential takes the 
form 



{7- 



where V{f3^) V{<p). Note that if the potential IS given 
by y = cA^ and if 7 = 4/3 (i.e., radiation), then the 
above Hamiltonian coincides with the one in general rela- 
tivity describing a massless scalar field and a cosmoiogicai 
constant A = nc. Thus intrinsically isotropic models (i.e., 
the isotropic models and the type I and V models, collec- 
tively characterized by V* = —6k when choosing a = 1 
for the type V models) are solvable if X(^nm) is chosen 
to depend only on since this leads to a 1-dimensional 
generalized Friedmann problem. Intrinsically isotropic 
stiff fluid models with the same scalar fleld potential, 
V = cA^, are also solvable since they lead to a separa- 
ble potential. The corresponding Jacobi metric therefore 
admits a second rank Killing tensor. 

a. The conformally coupled case. Conformal coupling 
correspond s to the choice A{(j)) = 1 — 0^/6, -6(0) — 1. 
Using Eq. (7.(;) leads to a redefinition of the scalar field 
given by the relation tp = v^tanh/S^ implying A (/?''') = 
cosh" ^ /?'''. For a model with a cosmoiogicai constant, a 
mass term, and a quartic term the total potential can be 
written in terms of the redefined fields as 

U = 24x(;l)[iy*e4^" -t- e«^"(Acosh4^t 

-H3kto^ cosh^ sinh^ + 3k\ sinh'' )] . (7.9) 

Obviously this is a SE-Hamiltonian. Apart from the 
"trivially" solvable intrinsically isotropic case with m = 
0, A = 0, A = (correponding to setting c = in the 
previous general discussion), there are also a number of 
other sets of values of the parameters for which the model 
is solvable [112|. The isotropic case with A = 0, A = 



was shown to be chaotic by Calzetta and El Hasi [113 



b. A solvable case with nonconformal quadratic cou- 
pling. Consider a nonminimally quadratically coupled 
model with A{(f>) = 1 — For these models it follows 
that the ^(0) = cA^ = k^^AA"^ case, discussed above, 
corresponds to a model with an arbitrary cosmoiogicai 



constant A, mass m = 2-\/— k^^A^, and a quartic term 
with A = 24k~^A^^. To have a physically reasonable 
mass term we must have A^ < 0. 

c. Brans-Dicke models. In this case the scalar field 
couphng is defined by the relations A{(j)) = ncj) , -8(0) = 
a;/0, and V{(f)) = 0. Redefining the scalar field hy (j) — 
= K-ii(/3t) where v = ^^3/(3-^ 2a;) leads to 
the Hamiltonian 



H 



+2'ix 



(nm) 



2-^/3^2) (7.10) 

3(2-7)/3»-(4-37)^'/3t" 



Note that this is a SE-Hamiltonian and that for radia- 
tion (7 = 4/3), it is equivalent to the general relativistic 
Hamiltonian with two noninteracting perfect fiuids, one 
stiff and the other radiation (provided one chooses a;(„,„) 
to be independent of (3^). Furthermore, if the fiuid term 
in the above Hamiltonian is zero then there is a 1-1 cor- 
repondence between solvable stiff fluid models in general 
relativity and vacuum solutions in Brans-Dicke theory. 

As seen from the above discussion there is a close math- 
ematical relationship between the nonminimally coupled 
scalar field Hamiltonians (and particularly Brans-Dicke 
theory) and the SE-Hamiltonians occuring in general rel- 
ativity. This explains the numerous exact solutions one 
has obtained in these theories and the equally numerous 
number of articles describing them in the literature (for 
Brans-Dicke theory see e.g., [114,115|). Moreover, the 
above discussion shows how one easily can find new ones 
if one is so inclined. 



C. A note on quantum cosmology 

The results of this article may be used as a first step in 
quantizing SH models. Ashtekar et. al. [116| have quan- 



tized the intrinsically multiply transitive vacuum models. 

and 



As can be seen from sections VI A 1 



VIA2 



these all 

have nonnuU decoupling in the Taub slicing gauge. For 
any nonnuU decoupling case the Hamiltonian takes the 
form 



iJ = ^i/^(y^p^) 



(7.11) 



when expressed in symmetry adapted dependent vari- 
ables an d ill a slicing gauge leading to decoupling. As 
done in ]116| ], one can make a canonical transformation 
such that each decoupled Hamiltonian is reduced to the 
square of a momentum P^, leading to 



(7.12) 



The only trace of the original potential is to be found 
in the ranges of the values of the new variables. At this 
stage one has a complete set of observables (constants of 
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the motion) a nd o ne can foUow the quantization proce- 
dure used in [|116| to quantize these models. However, 
even for nonnuU solvable models one obtains a complete 
set of observable s, an d the quantization procedure dis- 
cussed by Torre [1_17| can be used to quantize them. 



VIII. CONCLUDING REMARKS 

There are other methods than the ones presented in 
this article which exist for producing exact solutions. 
Their relationship to the present ones is discussed be- 
low, but the relationship is not completely understood 
and deserves further attention. The section is concluded 
with a general discussion on a number of different issues. 



constant of the motion linear in the first derivatives [120|. 
For Hamiltonian systems of this type it therefore seems 
reasonable that this symmetry corresponds to a Killing 
vector symmetry since the latter is associated with a con- 
stant of the motion which is linear in the momenta. How- 
ever, their method is also applicable to non-Hamiltonian 
problems. 

They apply their analysis to diagonal SH models. They 
rederive the Bianchi type I solutions with either a cos- 
mological constant or a perfect fluid and the orthogo- 
nal Bianchi type II stiff perfect fluid solutions. These 
Hamiltonian models do indeed admit Killing vector sym- 
metries. However, they also apply their method to a 
non-Hamiltonian tilted Bianchi type V stiff model an d 
thereby obtain an exact perfect fluid solution |100 102]. 



A. Relationship to other solution techniques 



3. Comparison with Hewitt's exact solution method 



1. Comparison with solution generating techniques 

Various solution generating techniques have been de- 
veloped for vacuum, electromagnetic or stiff perfect fluid 
spacetimes with one or two commuting Killing vectors 
(see e.g., [| ll^Hil )- These techniques rely on the exis- 
tence of symmetries which allow one to find new solutions 
from a given solution within the infinite-dimensional 
space of solutions being considered. For the finite- 
dimensional Hamiltonian problems studied here, one can 
also use the Killing tensor symmetries to generate new 
solutions from a particular one, but in practice this is a 
mute point since one finds the entire family at once. 

There are also solution generating methods which pro- 
duce new solutions from a particular one butwith differ- 
ent source or symmetry characteristics ji, IOC. 103]. Al- 
though the present approach analyzes separate Hamilto- 
nian problems, one could also use the variation of pa- 
rameters idea of ]103] to establish relationships between 
different Hamiltonian problems. 

For models with an infinite number of degrees of free- 
dom, one can impose conditions on various geometric 
quantities or on the functional form of the line element 
and still obtain a nontrivial problem corresponding to an 
invariant submanifold. This is in stark contrast to the 
situation for the finite number of degrees of freedom of 
the HH models where such conditions usually result in 
inconsistencies (except in the case when one has an un- 
specified function, like an arbitrary scalar field potential) . 
Thus it is critical to have systematic methods for finding 
invariant submanifolds for such models. 



2. Comparison with the exact solution method of Maartens 
and Wolfaardt 

Maartens and Wolfhaardt consider a certain class of 
systems of second order differential equations and find a 



For polynomial systems of ordinary differential equa- 
tions one can search for algebraic invariant curves, which 
then lead to exact solutions. Hewitt has applied such 
a method to 2-dimensional systems arising from the 
Einstein field equations for certain cosmological models 
]121], looking for linear and quadratic algebraic invari- 
ant curves. For 2-dimensional systems the existence of 
a sufficient number of such curves not only produces the 
corresponding exact solutions but also makes it possible 
to solve the full system. The search for these invari- 
ant curves is quite complicated and relies on algebraic 
computing, making it difficult to extend the approach to 
higher dimensions or to invariant curves of higher degree. 
Another consideration is the degree of the polynomials 
occurring in the system of equations, which must be suf- 
ficiently low for practical use. 

Such 2-diniensional polynomial systems can be derived 
if, for example, the problem is 2-dimensional and the 
Taub potential has at most two exponential terms. All 
the models of th is typ e which admit Killing tensors are 
given in section 



IV H 



The Killing tensors give rise to 
constants of the motion which are linear or quadratic in 
the momenta and involve exponential factors. 

Sometimes these constants of the motion lead to linear 
or quadratic curves in the polynomial system for certain 
values of those constants, but not always. Thus there is 
an overlap with Hewitt's method but it is not clear how 
large it is. So far all cases which have been found by 
that method correspond to the existence of Killing ten- 
sors although there are many Killing tensor cases which 
don't lead to linear or quadratic algebraic curves. On 
the other hand there may exist Hewitt cases which do 
not correspond to Killing tensors. 

The present approach has the advantage that one can 
immediately see whether or not exact solutions occur by 
inspection of a single function, the Taub potential, by 
hand, without attacking a whole system of differential 
equations using algebraic computing. Furthermore it is 
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easier to obtain an SE-Hamiltonian in '"standard" form 
than by rewriting the field equations in polynomial form. 
It also has the advantage that higher dimensional prob- 
lems and those with more complicated Taub potentials 
are easily handled. On the other hand Hewitt's method 
can be applied to 2-dimensional problems which do not 
arise from Hamiltonian systems of the usual type. 



Ha miltonian form in terms of the Misner parametrization 
[^,125, 126 1 and generalized as much as possible to the 
general case for all Bianchi types in Q. One can also 
use adapted variables to a ttempt a so-called regulariza- 
tion of the field equations |127.128]. Thus it seems clear 
that the tools presented here may prove useful in many 
applications involving the rich dynamics of HH models. 



B. Discussion 



The common practice of obtaining exact HH solutions 
in gravitational physics is to examine each new scenario 
as a new problem in isolation without considering its 
mathematical relationship to other such problems. How- 
ever, if a problem admits a Hamiltonian formulation, 
where the kinetic part of the Hamiltonian can be put in 
the "conformally flat" form, then the present approach 
can be applied. Since this approach has made it possible 
to unify, extend and bring order to many apparently un- 
related special results of the existing literature, it should 
prove to be a useful tool in future studies. 

There are many models not explicitly treated in the 
present article which could be investigated with these 
methods. Among these are a variety of static [Q and 
self-similar models [122 . Some timelike self-similar mod- 



els have already been treated in this way and some new 
solutions found |122]. Other sources or combinations of 
sources may also be considered leading to an abundance 
of models. 

All of these examples lie within conventional general 
relativity. However, the most likely applications will arise 
in exploring alternative gravitational theories. For exam- 
ple, of the numerous articles which regularly appear in 
this area, a randomly chosen one [124| analyzes a Bianchi 
type I supergravity model, which can be completely ex- 
plained in terms of the present analysis. This is not un- 
typical. We are not aware of any solvable case in the 
literature on HH models which cannot be explained by 
the existence of rank two Killing tensor and Killing vec- 
tor symmetries. It would be interesting to find an explicit 
solvable case solution not admitting such Killing symme- 
tries. 

Finally the present framework is not just valuable for 
the goal of searching for exact solutions but may serve as 
the starting point for a qualitative analysis of the more 
general behavior of the field equations. There are various 
kinds of Hamiltonian symmetries which may not be suffi- 
cient to lead to exact solutions. Nevertheless by adapting 
the variables to these symmetries, one obtains a simpler 
qualitative description. An example of such a symmetry 
is the homothetic Killing vector symmetry which many 
models exhibit ||2^]. This was exploited to develop an 
intuitive qualitative picture of the dynamics of a num- 
ber of models in This follows in the footsteps of 
the well known picture of the Mixmaster dynamics of 
the diagonal Bianchi type IX models when expressed in 



[8 

[9: 

[lo; 

[11 

[12; 

[is: 

[14 

[is: 

[16: 
[17: 
[is: 

[19 

[2o: 

[21 

[22: 



Supported by the Swedish Natural Science Research 
Council. 

Supported by the EC Human Capital and Mobility Pro- 
gram. 

C. W. Misner and A.H. Taub, Sov. Phys. JETP 28, 122 
(1969). 

M.P. Ryan, Jr. and L.C. Shepley, Homogeneous Rela- 
tivistic Cosmology, Princeton University Press, Princeton, 
1975. 

V.A. Belinsky, I.M. Khalatnikov, and E.M. Lifshitz, Adv. 
Phys. 31, 639 (1983). 

D. Kramer, H. Stephani, M. MacCallum, and E. Herlt, 
Exact Solutions of Einstein's Field Equations (Cambridge 
University Press, Cambridge, 1980). 

E. Verdaguer, Phys. Rep. 229, 1 (1993). 
A.H. Taub, Ann. Math. 53, 472 (1951). 

R.T. Jantzen, in Proc. Int. Sch. Phys. "E. Fermi" Course 
LXXXVI on "Gamov Cosmology" , edited by R. Rufiini 
and F. Melchiorri (North Holland, Amsterdam, 1987), 
p. 61, and in Cosmology of the Early Universe edited by 
R. RufRni, L.Z. Fang (World Scientific, Singapore, 1984), 
p. 233. 

C.W. Misner, K.S. Thorne, and J. A. Wheeler, Gravitation 
(W.H. Freeman and Company, New York, 1973). 
A.H. Taub and M.A.H. MacCallum, Commun. Math. 
Phys. 25, 173 (1972). 

A. Ashtekar and J. Samuels, Class. Quantum Grav. 8, 
2191 (1991). 

Y. Fujiwara, H. Ishihara and H. Kodama, Class. Quantum 
Grav. 10, 859 (1993). 

C.W. Misner, Phys. Rev. 186, 319 (1969). 

R.T. Jantzen Phys. Lett. 186B, 290 (1987). 

K. Rosquist, C. Uggla, and R.T. Jantzen, Class. Quantum 

Grav. 7, 611 (1990). 

R.M. Wald, General Relativity (University of Chicago 

Press, Chicago, 1984). 

K. Rosquist, preprint (1992). 

R.T. Jantzen, Ann. Phys. (N.Y.) 145, 378 (1983). 

M.P. Ryan, Jr., Hamiltonian Cosmology, (Springer, 

Berlin, 1972). 

B DeWitt, Phys. Rev. 160, 1113 (1967). 

C.W. Misner, in Magic Without Magic, edited by J.R. 

Klauder (Freeman, San Francisco, 1972). 

R.T. Jantzen and C. Uggla, Gen. Relativ. Grav. 24, 59 

(1992). 

A. Krasinski, Physics in an Inhomogeneous Universe, 
preprint (University of Cape Town, Cape Town, 1993). 



33 



To appear in Phys. Rev. D 



Preprint, March 1995, |gr-qc/9503061 



R.T. Jantzen, Phys. Rev. 37, 3472 (1988). 

K. Fiser, K. Rosquist, and C. Uggla, Gen. Relativ. Grav. 

24, 679 (1992). 

C. Lanczos, The Variational Principles of Mechanics, 
Fourth Edition (University of Toronto Press, Toronto, 
1970). 

R. Abraham and J.E. Marsden, Foundations of Mechan- 
ics, Second Edition, (Benjamin/Cummings Publishing 
Company, Reading, MA, 1978). 

C. Uggla, K. Rosquist and R.T. Jantzen, Phys. Rev. D 
42, 404 (1990). 

P.J. Giver, Applications of Lie Groups to Differential 

Equations (Springer, New York, 1986). 

K. Rosquist, J. Math. Phys. 30, 2319 (1989). 

W. Dietz, J. Phys. A. 9, 519 (1976). 

K. Rosquist and C. Uggla, J. Math. Phys. 32, 3412 (1991). 
H. Stephani, General Relativity (Cambridge University 
Press, Cambridge, 1990). 

C. Uggla, Class. Quantum Grav. 7, L171 (1991). 

C.G. Torre and I.M. Anderson, Phys. Rev. Lett. 70, 3535 

(1993). 

J. Wainwright, Gen. Relativ. Grav. 16, 657 (1984). 
R.T. Jantzen and K. Rosquist, Class. Quantum Grav. 3, 
281 (1986). 

C. Uggla and K. Rosquist, Class. Quantum Grav. 7, L279 
(1991). 

R.T. Jantzen, C. Uggla and K. Rosquist, Gen. Relativ. 
Grav. 25, 409 (1993). 

E. Kasner Amer. J. Math. 43, 217 (1921). 

E.R. Harrison, Mon. Not. Roy Astr. Soc. 137, 69 (1967). 

J. P. Vajk, J. Math. Phys. 10, 1145 (1969). 

J.L. Anderson, Principles of Relativity Physics (Academic 

Press, New York, 1967). 

W. De Sitter, Proc. Kon. Ned. Acad. Wet. 20, 229 (1917). 
P.T. Saunders, Non-isotropic Model Universes, Ph.D. the- 
sis, London, 1967. 

K.C. Jacobs, Astrophys. J. 153, 661 (1968). 
L Robinson, J. Math. Phys. 2, 290 (1961). 
A. Raychadhuri, Proc. Phys. Soc. Lond. 72, 263 (1958). 
A.G. Doroshkevich, Astrofiz. 1, 255 (1965). 
J.M. Stewart and G.F.R. Ellis, J. Math. Phys. 9, 
(1968). 

J. P. Vajk and P.G. Ehgroth, J. Math. Phys. 11, 
(1970). 

L. Hsu and J. Wainwright, Class. Quantum Grav. 3, 1105 
(1986). 

V. Joseph, Proc. Camb. Phil. Soc. 62, 87 (1966). 
G.F.R. Ellis and M.A.H. MacCallum, Commun. Math. 
Phys. 12, 108 (1969). 

A. Einstein, Preuss. Acad. Wiss. Berlin, Sitzber., 142 
(1917). 

J.J. Halliwell, Phys. Lett. 185B, 341 (1987). 
R. De Ritis, G. Marmo, G. Platania, C. Rubano, P. Scud- 
ellaro and C. Stornaiolo, Phys. Rev. 42D, 1091 (1990). 
R. De Ritis, G. Marmo, G. Platania, C. Rubano, P. Scud- 
cUaro and C. Stornaiolo, Phys. Lett. 149A, 79 (1990). 
G.F.R. Ellis and M.S. Madsen, Class. Quantum Grav. 8, 
667 (1991). 

G. Goetz, J. Math. Phys. 31, 2683 (1990). 

G. Goetz and D. Noetzold, preprint, Fermilab-Pub- 



1072 



2212 



[61 
[62 
[63 

[64' 
[65; 

[66 
[67 

[68' 
[69 
[70 
[71 
[72' 
[73; 

[74' 
[75; 

[76' 
[77 

[78' 
[79' 
[80 

[81 

[82' 
[83' 
[84' 
[85' 
[86' 
[87 
[88' 
[89; 

[90 
[91 
[92 
[93 

[94 

[95 

[96 
[97 



[99 
[100 

[101 

[102 
[103 
[104; 



89/235- A (1989). 

D. Noetzold, preprint, Fermilab-Pub-90/35-A (1990). 

D. R. Brill, Phys. Rev. 133B, 845 (1964). 

M. Cahen and L. Defrise, Commun. Math. Phys. 11, 56 
(1968). 

M.A.H. MacCallum, Gen. Relativ. Grav. 17, 659 (1985). 
M.A.H. MacCallum and S.T.C. Siklos, J. Geometry and 
Phys. 8, 221 (1992). 

G. Rosen, J. Math. Phys. 3, 313 (1962). 

J. Novotny and J. Horsky, Czech. J. Phys. B24, 718 

(1974). 

C.S. Kar, Phys. Z. 27, 208 (1926). 

G. C. McVittie, Proc. Roy. Soc. Lond. 124, 366 (1929). 
B.K. Datta, Nuovo Cim. 36B, 109 (1965). 

W.B. Bonnor, Proc. Phys. Soc. A 66, 145 (1953). 
K.C. Jacobs, Astrophys. J. 155, 379 (1969). 
V.A. Ruban, preprint 355 (A.F. loffe Institute, 
Leningrad, 1971). 

A. Barnes, J. Phys. A (1978) 11, 1303 (1978). 

R. Kantowski and R.K. Sachs, J. Math. Phys. 7, 443 
(1966). 

B. K. Datta, Ann. Phys. (N.Y.) 12, 295 (1961). 

K. Schwartzschild, Sitz. Preuss. Akad. Wiss., 424 (1916). 

F. Kottler, Annalen Physik 56, 410 (1918). 

H. Reissner, Annalen Physik 50, 106 (1916). 

G. Nordstrom, Proc. Kon. Ned. Acad. Wet. 20, 1238 
(1918). 

E. T. Newman, L.A. Tamburino and T. Unti, J. Math. 
Phys. 4, 915 (1963). 

R.C. Tolman, Phys. Rev. 55, 364 (1939). 

O. Klein, Ark. Mat. Astr. Fys. A 34, 1 (1947). 

G. M. VolkofT, Phys. Rev. 55, 413 (1939). 
M. Wyman, Phys. Rev. 75, 1930 (1948). 
K. Rosquist, in preparation (unpublished). 

H. A. Buchdahl, Astrophys. J. 140, 1512 (1964). 
H.A. Buchdahl, Astrophys. J. 147, 310 (1967). 

H.A. Buchdahl, Seventeen Simple Lectures on General 

Relativity Theory (Wiley, New York, 1981). 

W. Simon, preprint, (unpublished). 

A.B. Evans, J. Phys. A: Math. Gen. 10, 1303 (1977). 

C. B. Collins, Commun. Math. Phys. 23, 137 (1971). 
C.B. Collins and J.M. Stewart, Mon. Not. R. Ast. Soc. 
153, 419 (1971). 

C.B. Colhns, E.N. Glass, and D.A. Wilkinson, Gen. Rel- 
ativ. Grav. 12, 805 (1980). 

A.S. Kompaneets and A.S Chernov, Zh. Eksper. Teor. Fiz. 
47, 1939 (1964). 

M.A.H. MacCallum, Commun. Math. Phys. 20, 57 (1971). 

C. Uggla, R.T. Jantzen, K. Rosquist, and H. von Zur- 

Miihlen, Gen. Relativ. Grav. 23, 947 (1991). 

L Robinson and A. Trautman, Proc. R. Soc. A265, 463 

(1962). 

C.G. Hewitt, Class. Quantum Grav. 8, L109 (1991). 

J. Wainwright, W.C.W. Ince and B.J. Marshman, Gen. 

Relativ. Grav. 10, 259 (1979). 

M.A.H. MacCallum and A.H. Taub, Commun. Math. 
Phys. 25, 173 (1972). 

J. Wainwright, Gen. Relativ. Grav. 16, 657 (1984). 
R.T. Jantzen, Ann. Phys. (N.Y.) 127, 302 (1980). 
J. Barrow, Nature 272, 211 (1978). 



34 



To appear in Phys. Rev. D 



Preprint, March 1995, |gr-qc/9503d61 



[105] R. Maartens and S.D. Nel, Commun. Math. Phys. 59, 273 
(1978). 

D.L. Farnsworth, J. Math. Phys. 8, 2315 (1967). 
R. Maartens and M.F. Wolfaardt, Class. Quantum Grav. 
9, 1525 (1992). 

K. Rosquist, Phys. Lett. 97A, 145 (1983) (erratum in 
Phys. Lett. lOOA, 516 (1984)). 

K. Rosquist and R.T. Jantzen, Phys. Lett. 107A, 29 
(1985). 

R.T. Jantzen, Phys. Rev. D34, 424 (1986). 

D. N. Page, J. Math. Phys. 32, 3427 (1991). 
L. Bombelli and K. Rosquist, in preparation (unpub- 
lished) . 

E. Calzetta and C. El Hasi, Class. Quantum Grav. 10, 
1825 (1993). 

D. Lorenz, in Lecture Notes m Physics 205, edited by C. 
Hoenselaers and W. Dietz (Springer- Verlag, Berlin, 1984). 

E. Guzm'an, Astrophys. Space Sci. 173, 163 (1990). 
A. Ashtekar, R. Tate and C. Uggla, Int. J. Mod. Phys. 2, 
15 (1993). 

C.G. Torre, Class. Quantum Grav. 8, 1895 (1991). 
R. Geroch, J. Math. Phys. 13, 918 (1971). 
R. Geroch, J. Math. Phys. 14, 394 (1972). 
R. Maartens and M.F. Wolfaardt, Class. Quantum Grav. 
9, 1525 (1992). 

C. G. Hewitt, Gen. Relativ. Grav. 23, 1363 (1991). 

D. M. Eardley, Commun. Math. Phys. 37, 287 (1974). 
C. Uggla, Class. Quantum Grav. 9, 2287 (1992). 
L.O. Pimentel, Class. Quantum Grav. 9, 377 (1992). 
C.W. Misner, Phys. Rev. Lett. 22, 1074 (1969). 
C.W. Misner, in Relativity edited by M. Carmeli, S.I. Fick- 
ler, and L. Witten (Plenum Press, New York, 1970). 
K. Rosquist and R. Jantzen, Phys. Rep. 166, 89 (1988). 
C. Uggla and H. von Zur-Miihlen, Class. Quantum Grav. 
7, 1365 (1990). 



35 



To appear in Phys. Rev. D 



Preprint, March 1995, 



gr-qc/9503061 







I 


II 


VIo 


VIIo 


VIII 


IX 










1 


1 


1 


1 










-1 


1 


1 


1 







1 








-1 


1 





TABLE I. 

Canonical choices of the symmetry parameters n^°'\ For nondiagonal type II models the 
choice (ri(^\ n*^^-', n*^'^-') = (1,0,0) is more convenient than the one in the table. 



{ri,r2) 


A 


5 


(0,1) 


2q2 - qi 


q2 - qi 


(1,0) 


2qi - q2 


qi - q2 


(0,2) 




^{l2 - qi) 


(2,0) 


gi 


|(gi -92) 


(1,2) 




92 - qi 


(2,1) 


qi 


qi - q2 



TABLE II. 

The choices of (A, 5) for the two-term potential case 
which lead to a linear or quadratic potential. {ri,r2) 
are the new powers. 
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b 


S {P^,<lJ) 


(C,,dj) 


HKV KT case 


I 


6/1 


Pi = 2pj or 



Qi = 2qj 


Ci +di ^ 2{cj + dj) or 
Ci — di = 2{cj — dj) 


[A) 


II 




^ Pi =P2 or qi = q2 


ci - C2 = ±(di - d'z) 




111 


1 


1 Piq2 +P2qi = 


ciC2 — d\d2 =0 


(E) 


IV 




- Pi92 - P2qi = 


Cid2 — C2di — 




V 


1 

2 


Pi = 3pj and 

Z = ±1 

qj = Sqi 


Ci + di = 3(cj + dj) and 
Cj — dj = 3(ci — di) 


Z-[^ (C) 
(D) 


TABLE 111. 

Null and nonnuU parameters for 2-dimcnsional two exponential terms models admitting Killing tensors up to second 
rank. The index pair is a permutation of (1,2) where appropriate in the table. The parameters b and E are 
undefined in cases II and IV. The last column relates the different cases to the general HKV Killing tensor cases. 




Bianchi 
type 


dim 


vacuum A-term 


perfect 
fluid 


A-term plus 
perfect fluid 


I 


6 


Mink deS [§| 


FRW 


FRW-A 


I 


4 


Kasner Mink Saunders |@ 


Jacobs H 
Robinsoiil46| 
Raychaudhuri 1)471 
Doroshkevich [[fSf 
Ste-EU 01 
Vai-Eh ^ 


Saunders Q 


I 


3 


Kasner ||] Saunders @ 


Jacobs H 
Robinson[|46| 
Raychaudhuri p7[ 


Saunders Q 


V 


6 


Mink(Milne deS anti-deS 


FRW 


FRW-A 


V 


3 


Joseph H 


Ell-Mac (53I 




IX 


6 


— deS [|| 


FRW 


FRW-A (Ein Q) 



TABLE IV. 

SH models with 6-dimensional intrinsic symmetry group not including a scalar field. The abbreviations Mink, deS, 
anti-deS, Vai-Elt, Ste-EU, Ell-Mac and Ein stand respectively for Minkowski, de Sitter, anti-de Sitter, Vajk-Eltgroth, 
Stewart-Ellis, EUis-MacCallum and Einstein. The dimension column in this and all subsequent tables refers to the 
dimension of the spacetime symmetry group. 
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sym 
type 


dim 


vacuum 




A-tcrni 






em-term 


A-term -f 
em-term 


B I 


4 


Kasner M 

L—l' 


(e = ±l) 


Saunders Q 
N-H ^ie = 


1) 


-1) 


Rosen [p6| (e = —1) 
Kar ||§(e = 1) 
McVittie H (e = 1) 




B I 


3 


Kasner 


(. = ±1) 


Saunders Q 




-1) 


Datta ifife = -1) 
Bonnor ifrill (e = 1) 
Jacobs (e = —1) 




B II 


4 


Taub § (e 


= -1) 








Ruban |^ (e = -1) 
BcirnGS [|T4[] ^6 = ~1) 




B II 


3 


Taub § (e 


= -1) 












B III 


4 


K-S |7| {e 


= -1) 








Datta |||| (e = 1) 




KS 


4 


K-S |7| {e 


= -1) 












sss 


4 


Schwar. ^ 


1 (^=1) 


Kottler @ (e 






Reissner Q (e = 1) 


Nord. (e = 1) 


B VIII 


4 


Taub § (e 


= -1) 












B IX 


4 


Taub [i (e 
NUT jsil ( 


= -1) 


Brill |62| (e = 


-1) 









TABLE V. 

SH models with 4-dimensional intrinsic symmetry group not including a scalar field. The abbreviations B, SSS, 
N-H, Schwar., Nord. and NUT stand respectively for Bianchi, static spherically symmetric, Novotny-Horsky, 
Schwartzschild, Nordstrom and Newman-Unti-Tamburino. 
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r 


s 




equation of state 


Killing tensor type 


1 


- 


-aY* + bY^ 


a{5p + pf = 262(3p + p) 


null 


2 




-aY'' + bY^ 




null 


- 


1/3 


-aY'^ + b 


3p + p^2b 


null 


- 


2/3 


-a + bY-'^ 


p = a 


null 


1 


1/3 


aY^ ~ 2bY^ + Ea 


(5p + p-4Ea)2 8fe2 
2Ea — 3p — p a 


H-J (E = 1) 
harm. (E = -1) 


2 


2/3 


(6F)-^ [-c+(l + Yf + c_(l - Yf] 


p = c+(i + F)^ + c_(i-y)5 


H-J 






{6Yy^ ?lie[co{l + lYf] 


p = 9m[co(l + jy)^l 


harm. 



TABLE VI. 

Perfect fluids corresponding to Killing tensor cases. In the flrst two coulumns a "-" means that the value of r or 
s is arbitrary. The parameters a, &, c± are real constants while cq is a complex constant. The shorthand notation 
Y = e'^ is also used. For the last case, {r — 2, s — 2/3), it is not possible to write down an equation of state in closed 
form. Instead the expression for p{P^) is given in the fourth column. The abbreviations "H-J" and "harm." stand for 
Hamilton- Jacobi and harmonic respectively. The Hamilton- Jacobi Killing tensor case with (r, s) = (1,1/3) reduces to 
a Killing vector case when a — b. 
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Rianchi 



type 



dim 



1 < 7 < 2 



solution 



EPL 



Ortii. 



7 



flat FRW 



yes 



yes 



relevant sections 



^iA4|m| 



1 

4/3 

7 



Robinson 
Dor. 01 
JacobsM, S-E 



no 
no 
no 



yes 
yes 
yes 



VIAl 



IVE2 



IVE2 



IVE2 



m 
m 
[II 



Robinson_i|| 
Jacobs |4ql 



no 


yes 




VIAl 




IVE2 


no 


yes 




VIA] 




IVE2 



II 



7 < 2 
2 



C-S |9| 
Collins 



yes 

no 



yes 
yes 



VIA2 



VIA2 



|nj| 



II 



7 
2 

10/7 < 7 < 2 



Collins 
Collins 
Hewitt 



[32| 
[3S| 



no 
no 

yes 



yes 
yes 
no 



VIA2 



VIA2 



VIB3, 



Vl, |IV E 1 



IVE2I , Hlj 
■ 'III 



III 



1,4/3 

5/3 

2 




no 


yes 




VIA2 




IVEl 


no 


yes 




VIA2 




IVEl 


no 


yes 




VIA2 




IVE2 



KS 



1 

4/3 
5/3 
2 




no 
no 

no 
no 



yes 
yes 
yes 
yes 



VIA2 



VIA2 



VIA2 



VIA2 





IVEl 




IVEl 




IVEl 




IVE2 



V 



7 



open FRW 



yes 



VIAl 



[II 



V 



7 



E-M 



yes 



VIAl 



[VE2| , p] 



Vlh 



2(2cg+ l)/3 

7 < 2 

2 

2 

6/5 

2(2±cg)/3 
2(4cg - l)/3 
10/9 
2 



Collins 
Collins 
E-M H 
Collins 
Uggla 
U-R 1 37 1 
U-R ||__^ 
Wainlioil (h ^ -1/9) 
Wain, et al jlOO|(ft = 0) 



{h = 0) 
(ft/O) 
{h = -3/16) 



no 
yes 
no 

no 
no 
no 
no 

yes 
no 



yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
no 



VIA5 



VIA5 



via; 



via; 



via; 



via; 



via; 



IVEl 

EH 



IVEl 



IVE2, p] 



IVEl l 
VIB2I , 3[|li 



VIB3, [VEl 



IVE2 



IVEl 



VIII 



Jantzen 103] 



yes 



VIA2| , |[VE2| , p] 



IX 



closed FRW 



yes 



VIAl, [II 



IX 



Barrow \104] 



yes 



VIA2| , p/E% ^ 



TABLE VII. 

Hamiltonian spatially homogeneous perfect fluid models. The abbreviations Orth., S-E, C-S, K-C, "Dor.", K-S, U- 
R, R-U, E-M and "Wain." stand for Orthogonal, Stewart-Ellis, Collins-Stewart, Kompanets-Chernov, Doroshkevich, 
Kantowski-Sachs, Rosquist-Uggla, Uggla- Rosquist, EUis-MacCallum and Wainwright. One can set q = 1 in the 
expression cq — q/ y^p + 3(? for type Vlh- Note that when cq = 1/2, then VI/i=- 
implies that the model is orthogonal while a no implies a tilted model. 



=111. A yes in the "Orth." column 
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Bianchi 
type 


dim 1 < 7 < 2 






solution 






EPL 


Orth. 


V 


4 1 
2 


Farnsworth |l06| 
Maartens and Nel 


jios 


1 


no 

no 


no 

no 


V 


3 2 


Maartens and Wolfaardt 


107 1 


no 


no 


VIo 


3 4/3, 1.0411 < 7 < 1.7169 


Rosquist 


108 


, Rosquist and Jantzen |109] 


yes, yes 


no, no 




3 2 




Wainwright et al | 


100 




no 


no 


Vllh 


3 2 
2 




Barrow | 
Wainwri; 


W4\ 

iht et al 1 


io^ 




no 
no 


yes 
no 



TABLE VIII. 

Non-Hamiltonian spatially homogeneous perfect fluid models. The numerical values 1.0411 and 1.7169 are numerical 
approximations calculated by Rosquist and Jantzen. 
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